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INTRODUCTION 

The  problem  of  the  flexure  of  a  beam  under  the  action  of  a  single  load  stands 
today  in  much  the  same  state  as  Saint-Venant  left  it  half  a  century  ago,  nearly 
all  of  the  more  recent  work  having  been  devoted  to  the  problem  of  flexure  with 
continuous  systems  of  loading.  There  has  not  even  been  any  notable  additions  to 
the  forms  of  the  cross-section  of  the  beam  discussed  by  Saint-Venant  and  the 
rectangular  and  elliptical  cross-sections  may  be  taken  as  the  only  cross-sections  of 
simple  mathematical  form  for  which  full  solutions  exist.  It  appears,  for  instance,  that, 
although  the  general  solution  of  the  problem  on  Saint-Venant's  assumptions  includes 
the  case  of  a  beam  having  a  uni-symrnetric  cross-section,  no  such  case  has  been 
discussed  analytically,  far  less  numerically.  Yet  the  case  of  general  asymmetry  in  the 
cross-section  possesses  peculiar  interest  and  has  recently  become  of  real  importance 
in  connection  with  the  blades  of  aeroplane  propellers.  With  the  high  speed  of 
rotation  to  which  the  aeroplane  propeller  attains  there  is  a  strong  flexing  force  in 
the  air  resistance  and,  as  we  shall  see,  the  flexure  of  a  beam  of  asymmetrical 
cross-section  produces  a  torsion  of  a  by  no  means  negligible  amount.  To  the 
designer  of  the  aeroplane  propeller  this  is  of  the  utmost  importance  for  even  a  slight 
change  in  the  '  angle  of  attack '  of  a  propeller-blade  might  seriously  alter  the  effi- 
ciency and  speed  of  the  machine.  It  may  well  be  that  more  potent  causes  exist 
for  the  torsion  of  a  propeller-blade  subjected  to  a  continuous  system  of  pressure 
and  not  merely  to  a  single  end-load,  but  to  gain  some  knowledge  of  the  nature 
of  the  torsion  produced  by  flexure  in  the  simpler  case  of  the  present  paper  must 
be  regarded  as  the  first  step  towards  the  investigation  of  the  more  difficult  problem. 
There  is  the  further  limitation  to  the  generality  of  the  discussion  of  the  following 
pages  that  we  shall  treat  only  of  cross-sections  which  have  an  axis  of  symmetry 
perpendicular  to  the  plane  of  flexure.  The  problem  of  complete  asymmetry  still 
awaits  investigation. 

The  torsion  spoken  of  above  which  arises  in  the  flexure  of  a  beam  of  asym- 
metrical cross-section  can  in  general  terms  be  seen  to  arise  from  two  contributory 
causes.  Firstly  there  is  the  distortion  of  the  cross-sections  of  the  beam  which 
Saint-Venant  showed  to  be  a  necessary  accompaniment  of  flexure.  With  an 
asymmetrical  cross-section  the  distortion,  it  is  clear,  will  be  likewise  asymmetrical. 
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If  then  we  fix  the  end-section  in  almost*  any  way  which  will  allow  the  distortion 
of  the  theory  to  take  place,  this  distortion  and  particularly  that  part  of  it  which 
accompanies  the  curvature,  called  by  Thomson  and  Tait  anticlastic,  will  produce  a 
rotation  of  the  beam  as  a  whole,  together  with  what  we  may  call  an  '  impure ' 
torsion.  The  latter  we  shall  term  the  torsion  due  to  anticlastic  curvature.  It 
may  be  noted  in  passing  that. even  in  the  case  of  a  beam  of  symmetrical  cross- 
section,  there  will  arise  a  rotation  of  this  nature  if  the  end  fixing  be  asymmetrical 
in  character. 

This  torting  effect  of  the  distortion  of  the  sections  is  however  only  one  portion 
of  the  total  torsion ;  with  it  must  be  combined  the  effect  due  to  the  distribution 
of  shearsf  in  every  cross-section  of  the  beam — a  distribution  which  is  necessarily 
asymmetrical  if  the  cross-section  be  asymmetrical.  It  is  the  main  object  of  this  paper 
to  investigate  this  torsion  due.  to  shear  with  more  particular  reference  to  the  case 
where  the  cross-section  is  a  circular  sector,  complete  or  curtate.  The  solution  of 
the  corresponding  torsion  problem  exists  for  this  cross-section,  and  with  numerical 
reductions  such  as  we  should  expect  in  work  done  by  Saint- Venant  I.  Arithmetical 
work  is  even  more  necessary  in  the  flexure  problem  and  no  small  amount  of  care 
has  been  devoted  in  this  paper  to  make  clear  the  analysis  by  numerical  calculation, 
whenever  that  was  possible.  In  addition,  more  perhaps  with  a  view  to  strengthening 
the  authors'  faith  in  the  theoretical  results  than  to  proving  them  to  the  complete 
satisfaction  of  others,  a  series  of  rough  experimental  tests  were  carried  out  with 
various  beams  of  circular  sector  cross-section.  In  the  course  of  the  discussion  notes 
will  be  found  of  these  observations.  They  are  given  for  what  they  are  worth  and 
make  no  pretence  to  be  the  results  of  accurate  laboratory  experiments.  More 
painstaking  and  less  hurried  practical  investigation  must  be  postponed  to  times  of 
greater  leisure. 

§  1.     Description  of  the  contents  of  this  Paper. 

Saint- Venant  in  his  original  Memoir  §  of  1856  made  use  of  the  semi-inverse  or 
'  mixed '  method  of  treatment  and  assumed  both  that  the  stresses  perpendicular  to  the 
longitudinal  axis  of  a  flexed  beam  across  planes  parallel  to  that  axis  were  everywhere 
zero ||  and  that  the  extension  at  any  point  was  proportional  to  the  distance  of  that 

*  An  unimportant  exceptional  case  will  be  pointed  out  later  on,  §  9. 

t  We  shall  use  '  shear '  in  the  sense  of  a  stress  in  which  it  was  originally  introduced  into  elastical 
nomenclature  by  George  Stephenson,  the  corresponding  strain  being  called  a  '  slide '  by  Pearson.  Cf. 
Todhunter  and  Pearson's  History,  Vol.  I,  p.  882. 

\  Complex  Rendus,  T.  LXXXVII,  1878,  pp.  849 — 54  and  893 — 9.  The  paper  is  described  in  Todhunter 
and  Pearson's  History,  Vol.  n,  Part  I,  p.  191  et  seq.  Compare  also  Greenhill,  Mess,  of  Maths.,  Vol.  10 
(1880),  p.  83. 

§  Memoire  sur  la  flexion  des  prismes,  Journal  de  TAouville,  Tome  ler  (2e  seYie),  pp.  89 — 189. 

||  Or,  in  other  words,  that  there  was  no  mutual  stress  between  the  longitudinal  'fibres'  of  the  beam 
perpendicular  to  those  fibres.  This  is  the  phrase  used  by  most  writers,  Saint- Venant  among  them,  but  the 
suggestion  of  fibrous  material  which  it  conveys  makes  it  unsuitable  for  use  without  qualification  although 
it  brings  the  physical  meaning  of  the  hypothesis  well  before  the  imagination. 
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point  from  the  neutral  plane.  Clebsch*  showed  that  the  second  assumption  flowed 
from  the  first,  but,  as  his  equations  are  all  for  the  isotropic  case,  we  have  briefly  in 
§§  2 — 4  given  the  required  equations  generalised  for  the  aeolotropic  case  which  we 
treat  later  on  in  the  paper. 

In  the  second  part  of  the  paper  (§§  5 — 8)  the  general  equations — isotropy 
being  for  the  present  assumed — are  transformed  for  use  with  cylindrical  coordinates 
and  the  complete  analytical  solution  given  for  the  case  of  a  uniform  beam  with 
cross-section  bounded  by  two  concentric  arcs  and  two  radii,  the  plane  of  flexure 
being  taken  at  right  angles  to  the  radial  plane  of  symmetry  of  the  beam.  This 
plane  of  symmetry  is  therefore  the  '  neutral '  plane  as  regards  extension. 

Then  follows  (§§  9 — 14)  an  arithmetical  reduction  of  the  analysis,  especially 
in  connection  with  the  torsion  due  to  shearing  stress  and  its  dependence  on  the 
angle  and  on  the  curtateness  of  the  sector.  The  results,  both  of  theory  and  of 
experiment,  lead  to  investigation  of  the  distribution  of  vertical  shearing  stress 
along  the  central  axis  of  symmetry  of  any  cross-section  (§  12),  and  of  the  nature 
of  the  end-fixing  which  is  permissible  with  the  solution  obtained  in  accordance 
with  Saint- Venant's  assumption. 

So  far  the  discussion  in  connection  with  the  particular  case  of  the  circular 
sector  cross-section  has  been  made  on  the  assumption  of  isotropy  (and  of  uni-constant 
isotropy  so  far  as  the  numerical  work  was  concerned),  but  in  §§  15 — 19  there  is 
given  a  protective  method  which  leads  to  a  series  of  aeolotropic  solutions.  The 
numerical  work  based  on  this  brings  out  the  important  fact  that  aeolotropy  such 
as  occurs  in  a  wooden  beam  (with  the  grain  in  the  longitudinal  direction)  causes 
a  great  increase  in  the  amount  of  the  shearing  torsion.  It  is  found  that  the  torsion 
of  the  similar  isotropic  beam  has,  very  roughly,  to  be  multiplied  by  f-Ey/u.,  where  E  is 
the  longitudinal  stretch-modulus  and  p.  a  mean  value  of  the  corresponding  transverse 
slide-modulus.  For  very  high  values  of  E/p.  this  rule  ceases,  however,  to  be  even 
roughly  true. 

In  the  last  part  of  the  paper  will  be  found  brief  discussions  of  some  other  cases 
of  cross-sections  which  have  been  partially  worked  out.  To  begin  with,  an  investiga- 
tion is  made  of  the  algebraic  solution  of  Saint- Venant.  It  is  pointed  out  that  he 
appears  to  have  overlooked  the  torsional  effect  altogether  in  his  discussion  of  the 
forms  of  uni-symmetrical  cross-sections  which  arise  from  his  solution.  An  attempt — 
with  no  great  success — was  made  to  extend  his  method  to  cases  of  asymmetry. 
Then  follows  a  short  description  of  a  more  general  method  of  attacking  the  problem 
by  aid  of  suitable  conjugate  functions.  The  case  of  the  lemniscate  is  fully  worked  out. 
If  there  is  not,  as  we  fear  there  may  be,  some  slip  in  the  analysis,  this  section  presents 
the  interesting  feature  of  a  reversed  torsional  effect.  Some  details  are  given  of  a  few 
other  particular  cases.  We  hope  to  be  able  to  give  a  more  adequate  numerical  discus- 
sion of  these  cases  in  a  later  paper. 

*  Theorie  der  Elaslicitiit  fester  Kijrper,  S.  74.     Leipzig,  1862. 
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I.     GENERAL  THEORY  OF  THE  FLEXURE  PROBLEM. 

§  2.  General  expressions  for  the  displacements  in  an  aeolotropic  material 
wliere  there  are  no  mutual  stresses  between  the  'fibres '  perpendicular  to  those  fibres. 

Let  us  take  this  longitudinal  direction  to  be  the  axis  of  z  and  the  axes  of 
x  and  y  to  be  at  right  angles  to  it.  We  shall  denote  the  displacements  at  x,  y,  z 
by  u,  v,  iv  and  the  stresses  by  xx,  yy,  zz,  yz,  zx,  xy.  Assuming  uniform  aeolotropy, 
and  that  the  planes  of  the  coordinates  are  planes  of  symmetry  of  structure,  we 
shall  use  E  to  denote  the  stretch-modulus  in  the  direction  of  the  axis  of  z, 
/AJ,  /A2  to  denote  the  slide-moduli  in  the  planes  of  xz  and  yz  respectively,  and  77^  rj., 
the  Poisson's  ratios  for  contraction  in  the  directions  of  x  and  y  due  to  an  extension 
in  the  direction  of  z*. 

The  condition  that  the  stresses  perpendicular  to  the  longitudinal  axis  (z)  of  the 
beam  across  planes  parallel  to  that  axis  should  be  everywhere  zero  may  be  written 
analytically  as 

xx  =  yy  =  xy  =  0 
everywhere,  giving  immediately 

1  du          1  dv      dw' 
77,  dx          T;,,  dy      dz 


dx      dy 
The  body-stress  equations,  also,  reduce  to 

9  /du      dw\  _ 
dz  \dz       dx  I 


dw 


.(1). 


dz  \dz      dy 

8  /du      dw\  d  idv 

1  +        +^ 


.('2). 


Fixing  our  attention  on  w  to  begin  with,  we  deduce 

cfw 


•(3), 


dz3      dx*dz      dy'dz      dxdydz 

so    that   -  - ,    the   longitudinal    stretch,    is    linear   in    x,  y  and    in    z,   and  may  be 
written 


dw 


~  =  (A'x  +  B'y  +  C")  -  (Ax  +  By+C}z 


.(4). 


*  The  arrangement  of  the  subscripts  in  the  notation  is  made  different  from  the  usual  arrangement, 
because  so  far  as  the  aeolotropy  is  concerned  the  problem  is  essentially  two-dimensional. 
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Integrating  the  first  set  of  (l)  with  this  value  of  —  ,  we  get 


u  =  -  •*?!  CiM'a:2  +  ^y  +  C"«)  +  77!  (\Atf  +  Bxy  +  Cx)  z  +  F(y,  z) 

where  F,  G,  H  are  functions  to  be  determined  to  satisfy  the  remaining  equations 
of  (1)  and  (2). 

Finally  we  arrive  at  the  following  as  our  displacements, 
u  =  —  ryz  —  {•£-  fax'2  —  -r].2y')  A'  +  ^xyB'  +  rj^C'}  +  {^  fax"  —  r)ty")  A  +  ^xyB  +  ^xC}  z ' 

—  \A!z~  +  ^  Az3  —  yy  +  ftz  +  a' 
v  =  TXZ  —  {r),,xyA'  +  \  (f].2y"  —  vjiX?)  B'  +  r).2yC'}  +  faxyA  +  $  (r).,y~  —  f)iX"}  B  +  y^yO}  z 


B  C  (tf 

^  x*y  +  -  (-  + 

+  (A'x  +  B'y  +  C')  z  -  i  (Ax  +  By+C)  v  -  /3"x  +  a!'y  +  y' 

........................  (6), 

where  T  and  a,  /3,  y,  a!,  /8',  y',  a",  ft"  are  constants  of  integration,  and  X  ^s  a  function 
of  x  and  y  which  satisfies  the  equation 


<7>- 


These  equations  may  be  tested  by  substitution  in  the  '  equations  of  compati- 
bility.' They  are  a  generalisation  of  the  equations  given  'by  Clebsch*  for  the 
isotropic  case.  The  constants  a",  /3",  y  might  of  course  be  included  under  the 
arbitrary  function  ^.  but  it  is  clearer  from  the  physical  point  of  view  to  keep  them 
separate  since  the  linear  terms  of  the  displacements  are  those  which  must  be  used 
to  fix  the  end-conditions.  The  terms  in  y  give  rise  to  a  rotation  about  the  axis  of  z, 
and  the  other  notation  is  arranged  to  suit  the  other  possible  rotations  which  would 
be  given  by  a  =  a"  and  ft  =  ft". 

The  equations  (6)  and  (7)  give  the  most  general  displacements  possible  in  a 
material  which  may  be  supposed  to  consist  of  parallel  fibres  when  there,  is  no  mutual 
stresses  between  the  'fibres  '  in  any  direction  at  right  angles  to  the  fibrous  direction. 
This  was  the  assumption  first  made  by  Saint-  Venant  in  his  treatment  of  the  flexure 
of  a  uniform  prism  by  a  load  applied  at  its  free  end.  It  is  an  assumption  justified 
by  the  consideration  that  the  stresses  xx,  yy,  xy  must  all  be  zero  at  the  lateral  surfaces 
of  the  prism,  and  that  if  they  are  taken  to  be  zero  throughout  the  interior  of  the 
prism  we  can  deduce  a  solution  of  the  flexure  problem  which  can  be  made  to  satisfy 
certain  end-conditions  of  fixing  and  of  load.  It  is  true  that  in  this  matter  of  fixing 
and  loading  we  do  not  presume  that  the  distribution  of  stresses  over  the  cross- 

*  Theorie  de  I'Elasticite  :  Saint-Venant's  annotated  French  translation  (1883),  p.  150. 
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sections  at  the  ends  is  actually  attained  in  the  theory  as  the  experience  of  practice 
would  lead  us  to  desire,  and  that  we  have  to  make  use  of  the  principle  of  the  elastic 
equivalence  of  equipollent  systems  of  loads,  but  Saint-  Venant  on  physical  grounds, 
and  Boussinesq  by  analysis,  have  shown  so  conclusively  the  rapidity  with  which 
the  local  effects  of  external  forces  diminish  with  distance  from  the  point  of  application 
that  in  any  beam  having  a  length  several  times  greater  than  its  breadth  dimensions  the 
principle  may  be  freely  applied. 

We  proceed  now  to  the  particular  case  of  the  prismatic  beam  under  flexure,  and 
after  that  to  the  particular  prismatic  beam  which  has  for  cross-section  the  sector  of 
a  circle  or  annulus. 

§  3.  Particular  case  of  the  foregoing  theory  for  the  cylindrical  beam  under 
flexure. 

The  appropriate  direction  for  the  axis  of  z  in  the  case  of  a  right  prism  will 
of  course  be  that  parallel  to  the  generating  line.  We  will  take  the  end  that  is 
fixed*  to  be  z  =  0,  and  if  the  length  of  the  beam  is  I  the  free  end  will  be  z  —  l. 
We  assume  a  load  to  act  on  this  free  end  in  the  direction  of  the  axis  of  x,  and 
such  that  it  is  statically  equivalent  to  a  force  W  acting  through  the  centroid  of 
the  end  cross-section. 

Since  the  load  acts  at  right  angles  to  the  longitudinal  axis,  the  stress  zz  will 
be  zero  all  over  the  surface  z  =  I,  that  is 


giving  A'  =  Al,        B'  =  Bl,         C'=Cl, 

and  therefore  d-^  =  (Ax  +  By  +  C)  (l-z)  ...........................  (8). 

' 


Now  let  us  consider  the  conditions  affecting  the  resultant  traction  over  any 
cross-section  of  the  beam.  These  conditions  are  obtained  by  forming  the  equations 
of  equilibrium  of  the  portion  of  the  beam  between  the  cross-section  in  question  and 
the  end  z  —  l,  the  only  forces  being  the  load  W  and  the  tractions  over  the  cross- 
section.  The  equations  are 

W  =     I  zxdxdy,     0  =  I  I  yzdxdy,     0=1    zzdxdy, 

—  Wy=\\(x.yz  —  y.xz}dxdy,      —  W(l  —  z)=  I  \x.  zzdxdy,     0  =  I  \y  .  zzdxdy 

.......................  ,(9), 

y  being  the  y-coordinate  of  the  centroid  of  the  cross-  section  and  the  integrations 
taking  place  over  the  area  of  the  cross-section. 

*  The  exact  mode  of  fixing  the  end  of  the  beam  and  of  applying  the  load  at  the  other  end  will  be 
discussed  later  (§  14). 


or 
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Consider  first  the  integrals  concerning  zz. 
From  0  =  I  I  zzdxdy   and  equation  (8)  we  derive 

0=11  (Ax  +  By+C)  dxdy, 

0-(t(A(x 
-))(A(x- 


where  (x,  y)  is  the  centroid  of  the  cross-section,  and  for  this  equation  to  be  satisfied 
we  must  have 

Ax  +  £y+C=Q    ..............................  (10). 


or 

we  have 


From  0  =     I  y  .  zzdxdy, 

0  =  I  I  (Axy  +  Bif  +  Cy)  dxdy, 


+  Axy  +  By2  +  Cy)  dxdy, 
and,  using  (10)  and  the  property  of  the  centroid, 

0  -      (A  (x-x)(y-y)  +  B(y-  y)*)  dxdy. 


If*  therefore  we  take  the  axes  of  x  and  y  to  be  parallel  to  the  principal  axes 
of  the  section  through  the  centroid  and  thus  make  the  product  moment 

\\(x-x)(y-y)dxdy 

vanish,  we  must  have 

B  =  Q     ....................................  (11). 

Again  from 

-  W(l  —  z)  =  I  Ix.  zzdxdy 

and  equation  (8)  we  have 

-W(l-z)=E(  ((Art  +  Bxy  +  Cx)  (I  -  z)  dxdy, 

and  using  (10)  and  (11) 

-W=E((A(x-x)tdxdy, 

W 
A=~EI 


where  /=  I  I  (x  —  xfdxdy, 

the    second    moment    of   the    cross-section    about    the    line    through    the    centroid 
perpendicular  to  the  plane  of  flexure. 

*  It  is  of  course  not  necessary  to  make  this  simplification. 
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With  these  simplifications  the  longitudinal  stretch  is  now  given  by 

dw          W,         .,,      . 
-  =  --(x-x)(l-z), 

or,  if  we  select  as  origin  some  point  on  that  principal  axis  through  the  centroid 
which  is  at  right  angles  to  the  plane  of  flexure,  so  that  x  =  0, 


An  examination  of  the  values  of  u,  v,  w  in  (6)  will  show  that  such  a  change 
of  origin  will  be  absorbed  by  the  arbitrary  constants  and  the  function  x- 

The  expressions  for  the  displacements  now  become 


W 

''  =  ~  Tyz  +  ~ 


w 

V  —  TZX  +  -J=TT 


...(14), 


where  x  (x>  y]  satisfies 

-0    ..............................  (15). 


To  apply  the  remaining  traction  conditions  of  (9)  we  require  expressions  for 
the  shears  xz  and  yz.     These  are 


W  (E—i)tiL\ 

~m 

(16). 


The  first  of  (9)  may  be  written 


w 

W=     \zx 


\         \dx       dy        I  )    ) 
from  (13)  and  the  third  body-stress  equation  of  (2), 

C  f  '  rt  -^-s  7]  x~s  \ 

or  W=W+\\(—(x.zx)  +  —  (x.yzY\dxdy 

j  j  \dx  ^  dy^          '  / 

—  W+  Ix  (xz  cos  (xv)  +  yz  cos  (yv)}  ds 

by  Green's  theorem,  v  being  the  normal  outwards,  the  integration  being  round  the 
boundary  of  the  cross-section,  and  in  order  that  this  may  be  true  we  must  have 

xz  cos  (xv)  +  yz  cos  (yv)  =  0  ...........................  (17) 

all  over  the  boundary.     We  shall  refer  to  this  as  the  Boundary  Condition. 
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Again  substituting  from  (16)  in 

r  r  ^ 

0  =  I  \yz  dxdy 
we  get  0=     fta*       ^      L£jfeft    xy-a  +  a"dxdy. 


The  integral  I  \xdxdy  vanishes  because  of  our  choice  of  origin,  and  I  \xydxdy  will 

vanish  *  if  the  cross-section  be  symmetrical  about  the  axis  of  y.  This  is  true  in 
all  the  cases  treated  in  this  paper,  and  we  will  make  this  additional  assumption 
in  all  that  follows.  We  thus  have 


and,  transforming  by  Green's  theorem  into  an  integral  round  the  boundary  of  the 
cross  -section,  we  get 


os(.^)cfe     .....................  (18), 

an  equation  which  we  shall  refer  to  as  Condition  I. 

From  the  last  remaining  of  the  traction  conditions  (9) 

f  [       "•"          •"> 
—  Wy  =  I  \(x  .  yz  —  y  .  xz)  dxdy, 


we  derive 


This  we  will  refer  to  as  Condition  II. 

We  have  now  reduced  the    physical  problem  to  the  mathematical  problem  of 
finding  the  function  x  (x>  y)  which  will  satisfy  the  differential  equation 


and  at  the  same  time  satisfy  the  three  conditions  given  in  (17),  (18),  (19)  for  the 
particular  cross-section  we  have  in  mind.  This  at  least  is  the  ideal  problem  ;  in 
reality  the  more  immediate  problem  is  to  find  the  particular  cross-section  which 
will  possess  as  nearly  as  possible  the  special  characteristics  we  want  to  investigate 
and  which  will  at  'the  same  time  permit  of  a  suitable  simple  function  x  being 
found  for  it. 

In   addition   to    our   original    assumption    of   no    mutual    stress    between    the 

*  The  integral  will  also  vanish  if  the  origin  be  taken  at  the  centroid  but  soluble  cases  will  then 
usually  be  symmetrical  about  both  principal  axes. 
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longitudinal  'fibres'  of  the  beam  perpendicular  to  their  direction,  it  is  to  be 
remembered  that  we  have  also  assumed 

1°,  that  the  plane  of  flexure  is  parallel  to  a  principal  plane  of  the  beam  through 
the  centroidal  line,  and 

2°,  that  the  beam  is  symmetrical  about  the  plane  through  the  centroidal  line 
and  perpendicular  to  the  plane  of  flexure,  an  assumption  which  makes  1°  imme- 
diately satisfied. 

§  4.      The  case  of  isotropy. 

It  will  be  convenient  at  this  stage  to  postpone  the  aeolotropic  discussion  and 
to  summarise  the  foregoing  for  the  isotropic  case  where  77,  =  7/2  =  17  and  ^  =  /A,  =  /*. 

We  have  taken  the  longitudinal  axis  of  the  beam  to  be  the  axis  of  z  and  the 
plane  yz  to  be  a  plane  of  symmetry,  the  flexion  taking  place  in  the  plane  zx.  The 
fixed  end  of  the  beam  is  at  z  =  0  and  the  loaded  end  at  -z  =  I. 

The  displacements*  are 

W 

- 


(20), 


W 
v  =  TZX  +  -j-rf^xy  (l  —  z)  +  yx  —  az  +  ft' 


W 
iv  =  x (x,  y)-[x (lz -  iz2)  +  xy*-\ - ft'fx  +  a"y 


and  T,  a,  ft,  y,  a!,  ft',  y',  a",  ft"  are  arbitrary  constants. 

In  addition  we  have  the  following  conditions  to  be  satisfied. 
(i)     The  Boundary  Condition  : 

xz  cos  (xv)  +  yz  cos  (yv)  =  0    ..................  .  .....  (22), 

•~-         (  r)v      W 

where  xz  =  p     -ry  +      -        (!^'  +  (l  ~  £*?)  ^)  +  £  ~  £" 


and 


(ii)     Condition  I  : 

0  =  |(X  -(«-«'»  ^s  (yv)d8    .....................  (23), 

the  integral  being  round  the  boundary  of  any  cross-section,  and 
(iii)     Condition  II  : 


W 

............  (24)- 


*  Compare  Love's  Elasticity,  2nd  Edition,  Chapter  xv,  Equations  (12),  p.  319.     The  axes  of  coordi- 
nates and  notation  used  in  the  present  paper  make  comparison  easy  with  the  results  given  in  Chapter  xv. 
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II.     CASE  OF  A  BEAM  HAVING  AS  CROSS-SECTION 
A  CIRCULAR  OR  ANNULAR  SECTOR. 

§  5.      The  equations  suitable  to  the  case  of  the  circular  sector  cross-section  with 
isotropy. 

For  our  problem  of  the  sector-shaped  cross-section  we  must  now  transform 
our  equations  to  cylindrical  coordinates 
z,  r,  6.  Let  0,  the  centre  of  the  circular 
arcs  be  our  origin,  and  let  the  axis  of  y 
be  the  line  6=0.  We  shall  take  the 
sector  to  be  bounded  by  the  arcs  r  =  a 
and  r  =  aa,  and  the  lines  Q—  +  Jy. 

The  boundary  condition 

xz  cos  (xv)  +  yz  cos  (yv)  =  0 

is,  on  substitution  of  xz  and  yz, 

9v        /     \     9v 
-  """  ' xv)  +  -—  cos  (yv) 


FIG.  1. 


W 
=  r(y  cos  (xv)  -  x  cos  (yv))  +  «r[(J 


i  cos  (xv)  +  (2  +  7))  xy  cos  (yv)] 


>;») (25), 

and  we  can  now  obtain  the  conditions   that   are  to  be    satisfied    on    each  of  the 
edges  of  the  sector. 

1°.     Over   the    external  arc,    where   r  =  a   and    —  ay^^  +  Jy,  the  boundary 
condition  becomes 

[2»]       —  r(acos  6  sin  6  —  a  sin  6  cos  6)  +  -prr[(^rta^sin-0  +  (l  —  AT?)  crcos2#)  sin  6 

\dr/r=a  &1 


or 


+  (2  +  TJ)  a2  sin  0  cos2  ff]  -  (/8  -  ft")  sin  6  +  (a  -  a")  cos  6, 

iv)     =  '~M  ^ + ^ sin  e + * sin  3^  "  ^  ~  ^ sin  * + (a  ~  a") cos  e (26)' 

2°.     Over  the  internal  arc,  it  is  easily  found  that 

(26)'. 


3°.      Over   the   radial    edge,    o0^r<a,    and    #  =  Jy,   the    condition    is,    in    the 
same  way, 

Wr2 

(*  ~  W cos  *r + * cos  W 

-(/3-  ft")  cos  ly  -  (a  -  a")  sin  Jy (27), 
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and,  4°,  over  the  other  radial  edge  0=  —  ^y> 

=  TT  +     [(i  ~  |7?)  cos  *r  +  f  cos  f  r] 


(a-a")siniy  ......  (27)', 

the  only  difference  between  3°  and  4°  being  the  sign  of  the  term  containing  a  —  a". 
But  if  x  is  to  be  one-valued,  the  conditions  3°  and  4°  cannot  both  be  satisfied  unless 
a  =  a".  It  will  be  evident  after  we  have  found  x  that  the  physical  interpretation 
of  this  is  that,  in  a  beam  which  is  symmetrical  about  the  plane  of  yz,  there  is  no 
lateral  slide  when  a  flexing  load  is  applied  in  the  principal  plane  zx. 

§  6.     The  solution  for  x  suitable  to  the  case  of  the  sector. 
We  shall  now  assume*,  as  our  solution  of 

^  +  ^_ 

9x2  •    Sf  ' 


n30   ...(28), 

and  proceed  to  determine  alt  a,,  a3,  A,,  A,,  ...,  Bl}  B^,  ...  and  the  appropriate  values 
of  m  from  the  conditions  obtained  above. 

From  the  boundary  condition  over  the  straight  edges 

a    \  /  7?    \ 

3L  =  2  (  AnT™'1  +  -—-.  }  m  cos  Amy  +  a,  cos  \y  +  a.2r.2  cos  y  +  a3r~  .  3  cos  f  y 

f      I 

y    ......  (29), 


which  will  be  satisfied  only  if 

my 

2  = 

where  i  is  zero  or  a  positive  integer, 

T 
and          a^-tf 

Hence  we  have 


cos  — 


*  We  have  to  make  x  an  odd  function  of  6  to  satisfy  the  equations  of  §  5.  Physically  this  means  that 
the  longitudinal  displacements  in  the  beam  are  symmetrical,  with  change  of  sign,  about  the  axis  of 
symmetry  of  the  cross-section. 
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Applying  now  the  boundary  condition  over  the  arcs,  we  have 


r-a    i= 


I  <ain 

sin 


Wa?  \~  1 

-gT    (f  +  ^)sin^  +  f  sin 30    -(£-£") sim 

**•*    L  J 


or 


along  with  a  similar  equation  in  «„  for  the  internal  surface. 

These  we  may  write  for  brevity  as 

t 

00 

2   Cm  sin  m0  =  2cp  sin  p^          O=l,2,  3]  ...............  (32), 

»=o  P 

where  cp  is  known  and  (7m  contains  the  quantities  Am  and  .Sm  which  are  to  be  found 

We  can  now  determine  the  constants  Cm  by  Fourier's  method  and  thereafter 
obtain  Am  and  Bm.  Multiply  both  sides  of  the  equation  (32)  by  sin  mO  and  integrate 
between  the  limits  —  \y  and  |y.  Then 


fiv  (\t 

m\       s\tfm.ddd  =  '2,cp\       smp 
f-tn  P     J  -47 

fi? 

=  ^cp  I        {cos  (p  -m}6  —  cos  (p  +  m)  6}  dO, 
p     }  -4v 


and 


.................................... 


y  p       m  —  _ 
Replacing  the  notation  Cm  and  cp,  we  have,  after  a  little  reduction, 


-  L 


"         l       j     "' 
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the  solution  of  which  is 

A 


(ma  -  4)  my    a™  -a 


_ 

m'2-l          m2-9 


D 


(m2  —  4)  my 


-(-I)'' 


co8 


3(1- 
m3-9 


a™  -a™ 


El    my     \m?- 

We  thus  have  x  and  therefore  the  displacements  u,  v,  w  expressed  in  terms  of 
and  the  arbitrary  constants  a,  /3,  /3",  y,  a',  /3',  y' : 


W 


w—\  — 


m 


8r  (  -  I)1'     a" 

~  ~ 


s  -  4)  my  ~cCM  -  of* 

W  cos       /3  +  27}      3(1- 


-  a 


my    \ma-l         m2-9 

T  V   ^  cos  sr  /3  +  ^ 
~-~~ 


g«»-gfl- 


•     mft 
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2  »***  **V    T    r»   77  J"  I     ^ 

cos  y  3ii/  L 

~-WT\-X  (lz  ~  K)  +  ^1  ~  F'x  +  a'V  +  / 


cos  f  y 


.(35), 


where 
and 


;  =  rsin0,    y=rcos0 


7T 

y' 


We  note  that  the  displacement  in  the  direction  of  the  flexing  force,  u,  is 
symmetrical  both  in  sign  and  magnitude  about  the  radius  of  symmetry,  while  v 
and  w  have  both  alternating  symmetry. 

§  7.  Application  of  the  Conditions  I  and  II  and  derivation  of  the  value  of 
the  torsion  T. 

We  have  now  to  apply  the  Conditions  I  and  II,  equations  (23)  and  (24)  of  §  4, 
and  will  thus  determine  two  relations  between  T  and  /3  —  /3"  from  which  these  may 
be  obtained. 
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(I)  Condition  I,  that  is 

fxcoa(yv)ds  =  0    ..............................  (36), 

the  integration  being  round  the  contour  of  any  cross-section. 

Taking  the  four  parts  AB,  BC,  CD,  DA  (see  Fig.   1)  separately,  we  have: 

1°.     Along  AB,  cos(yv)  =  —  sin^y,  and 

fs  ra 

X  cos  (H  ds  =  ~  sin  iy       Xe=iy  dr- 

J  A  J  Og 

2°.     Along  CD,  in  the  same  way, 

en  fao  rn 

X  cos  fa*}  ds  =  ~  sin  ir       Xe=  -  ^  (  ~  dr)  =  ~       X  cos  (yv)  ds. 

J   (-'  J  a  J  •"- 

3°.  Along  BC,  cos(yv)  =  cos0,  but  x  is  °dd  about  0  =  0,  and  therefore  the 
integral  vanishes. 

4°.     Along  DA,  the  integral  again  vanishes. 
Thus  Condition  I  is  identically  satisfied. 

It  is  easy  to  see  from  (36)  that  Condition  I  is  identically  satisfied  whenever  the 
section  has  uniaxial  symmetry  about  an  axis  perpendicular  to  the  plane  of  loading. 
For  in  this  case  x  changes  sign  as  we  cross  the  axis,  and,  if  P  and  P'  be  corresponding 
points  on  either  side  of  the  axis,  an  element  ds  at  P  will  correspond  to  an  element  ds 
at  P',  cos  (vy)  will  be  the  same  at  both  points,  but  x  will  change  sign  as  containing 
only  odd  powers  of  x,  and  thus  the  two  elements  of  the  integral  in  (36)  just  cancel. 

(II)  Condition  II,  that  is 


W  W 

-2l(2  +  kil)tfy  +  E 

When  the  integrations  are  carried  out,  we  obtain 


Tjr  4  a3  —  a03siniy      2W/J.,  xrn  .    . 

3a*-df     y      ~  5^7  ^      a")Wi-^)si 

arid,  transforming  the  surface  integral  by  Green's  theorem, 


Taking  the  line  integral  in  four  parts,  as  in  the  previous  case,  we  have  : 
1°.     Along  AB,       cos(yi>)=  —  sin^y,     cos  (xv]  =  cos  ^  y,  and 

tK  fa  fa 

x(xcos(yv)-ycos(xv))ds=  -       Xa=i-,  (r  sin2  2  7  +  r  cos2^y)cZr=  -       Xe= 

J  A  J  a0  J  a0 


3 
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2°.     Along  CD,  in  the  same  way, 

rn  fa 

I  x (x cos d/v) ~ y cos (xv)) d*>  =  +    x»=A  r(~ dr), 

Jo  J  a0 

and  therefore  the  sum  of  the  values  along  AB  and  CD 

_  o 

—       —i 

3°.     Along  BC, 

cc  r 

X  (x  cos  (yv)  —  y  cos  (xv))  ds  =  I  x(a  sm  ^  cos  ^  ~  a  cos  ^  sin  ] 
J5  J 

and.  4°,  along  DA  the  integral  again  vanishes. 

ra 
Thus  the  whole  integral,  P,  is  simply  —2  I    X0=A   rc^r< 

J«o 

Substituting  x  from  (30)  we  find  : 


P=  -2 


f 


73 


_  r» 


_  sn  -   niv 
-2  ^  U2    7 


-     sn 


tan 


4  -  a«4) 


^4m  and  .Bm  being  given  by  the  formulae  of  equation  (34),  and  on  the  substitution  of 
this  value  of  P  in  (37),  of  the  values  of  Am  and  Bm  from  (34)  and  the  value  of  /  we 
obtain,  remembering  that  E=  2/z.  (1  +17),  after  some  reductions  the  following  equation 
for  T  : 

(,v«t+2  .*,  m+2\2 

~ao        ) 


t00  fi4 

->4-«04)(tany-y)+2   /    »     >,\  /        ^Y 
i=0(m2-4)  (m  +  2) 


my 


64 


itoK-4)(m-2)my          «*"-<»„"" 

/4(3  +  27?)  _  12(l-2iy)\  (ctm+3 - a0m+8)  (am+2 - abm+i') 
V  m2-!  m2-9     J  "  a"" -a,"" 

V      2cos|y     /4(3  +  2T7)       12(1-27?)^  etX5(^m'3-<''3)(«m'2-<""2) 
m3-!  m2-9     /  a"" -a,"" 


+  T^  (1  -  2>?)  cos  ly  tan  f  y  (a5  -  a/) 


,  ...(39). 


The  value  of  /3  —  /3"  still  outstanding  will  flow  from  the  nature  of  the  terminal 
conditions. 

As  we  shall  see  in  §  9,  it  is  this  quantity  r  which  gives  the  torsion  in  the  beam 
due  to  the  asymmetrical  distribution  of  the  shears  in  any  cross-section. 
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§  8.     The  solution  in  the  case  of  the  complete  or  non-curtate  circular  sector. 

For  the  case  of  the  complete  sector,  that  is  when  a0  is  taken  to  be  zero,  the 
solution  of  §§  6  and  7  simplifies  into  the  following  : 


v  =  S^mrm  sin  m0-(/3-  /3")  r  sin  6  +  ~  —  r2  sin  20 

2  cos  y 


where      A  "1)'    J_  4.  (.^IL^  cosiy  /3  +  ^  _  3(1-2 

'"—^       El          my    \m?-l~      m2-9    /a 


1  T 

and  m  =  ( 2i  +  1 )  —  , 

7 
and  the  torsion  is  given  by 


r     — 


I      2  cos  ^y     /4  (3  +  2ij)  _  12(l-2q) 
El     .-.o  (m  +  2)my  \    m*-l  m2-9 


^y  tan  |y-f  (1  -2^)  sin  iy-f  (l+^)l  -  sin  | 

(41). 


III.     NUMERICAL  WORK  IN  CONNECTION  WITH  THE  TORSION  DUE  TO  FLEXURE. 
DISCUSSION  OF  SHEAR  DISTRIBUTION  AND  OF  TERMINAL  CONDITIONS. 

§  9.      Torsion  of  beam  with  Jlexure. 

As  stated  at  the  beginning  of  the  paper  the  torsion  of  the  beam  about  its  axis 
is  due  (i)  to  the  resultant  torsional  couple  which  arises  from  the  shears  in  directions  at 
right  angles  to  the  axis  and  (ii)  to  the  anticlastic  curvature  caused  by  the  inequality 
of  the  longitudinal  tension  and  consequent  inequality  of  the  Poisson  contraction  over 
any  cross-section.  For  each  of  these  effects  to  cause  torsion  it  is  of  course  essential 
that  the  cross-section  be  asymmetrical. 

To  determine  the  amount  of  torsion  let  us  take  the  vertical  displacements  of 
points  on  the  horizontal  axis  of  symmetry  at  the  fixed  end  and  at  the  other  end 
of  the  beam.  These  are  given  by  the  equation  (20) 

W 

-  z)  (x-  -  if]  +  ^  -  ±z*}  -yy  +  $z  +  a!. 


Firstly,  at  the  fixed  end  of  the  beam  the  displacement,  (i)  at  the  outermost  point 
of  the  axis  of  symmetry,  namely  (0,  a,  0),  is 

W 


3—2 
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and  (ii)  at  the  inner  end  of  the  axis  (0,  o0,  0), 

W 


W 
giving  ua  -  Ua0  =  -  jf-  ^r)l  (a2  -  a02)  -  y  (a  -  a0). 


The  first  term  arises  through  the  anticlastic  curvature  and  the  second,  being 
merely  a  rotational  term  affecting  the  whole  beam,  depends  on  the  particular  method 
of  fixing  the  end  of  the  beam.  When  we  consider  the  torsion  of  one  end  relative 
to  the  other,  it  will  cancel  out. 

At  the  end  where  the  load  is  applied  the  corresponding  displacements  are 

W 

(i)  at  (0,  a,  1)  ua'  =  -  -ral  +  -     ^  -ya  +  (3l  +  a', 


W 

and  (ii)  at  (0,  aa,  1)  uao'  =  -  raj  +        ±l3  -  yac  +  ftl  +  of, 


giving  Ua  -  Ua0'  =  -  rl  (a  -  a0)  -  y  ( a  -  a0), 

W 
the  terms  -^rr  \ Is  +  ftl  +  a!  disappearing  since  they  are  really  the  expression  for  the 


flexure. 

Taking  fc»  =  U  =  -          y  («  +  O  -  y 


—  Ua 


,  u,a   —  u,a  j 

and  v>i    =  -  —  —  Ti  —  y, 

a-a0 

as  approximately  giving  the  small  angles  of  rotation  at 
the  two  ends  of  the  beam,  we  have  the  total  torsion  - 
of  the  loaded  end  relative  to  the  fixed  end  about  the 
longitudinal  axis  of  the  beam 

W 

Direction  of 
the  direction  of  the  torsion  \\il  —  -^tr>  being  as  shown  in  Flexure 

the  diagram  of  the  free  end  of  the  beam  (Fig.  2).  Fl°-  2- 

In  connection  with  the  anticlastic  curvature  it  is  to  be  noted  that  it  is  zero 
at  the  free  end,  there  being  no  stretch  there,  and  that,  accordingly,  there  is  no 
distortion  of  the  surface  at  the  free  end  except  for  the  Saint- Venant  distortion  in 
the  direction  of  the  beam's  longitudinal  axis.  Nevertheless  if  we  were  to  clamp 

*  At  first  sight  this  expression  for  i/f0  looks  strange  in  that  it  would  appear  that  the  angle  of  torsion 
due  to  the  anticlastic  curvature  would  be  infinitely  large  for  the  limiting  case  of  a  rectangular  beam 
which  is  given  by  taking  the  curtate  sectors  to  have  infinite  radii  a  and  a0  but  such  that  a  —  «0  remains 

W 

finite.     In  taking  the  limit  of  -=,-=(«  +  aA  however,  it  is  to  be  remembered  that  we  must  take  the  limiting 

&  I 
expression  for  /  as  well,  and  7  being  \  (a*  -  a04)  (y  -  sin  y)  it  is  seen  that  the  angle  of  torsion  is  proportional 

to  ;-j jT-/ r  which  gives  zero  angle  of  torsion  for  the  rectangular  beam  as  it  should  do. 
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the  fixed  end  in  such  a  way  that  the  ends  of  the  axis  of  symmetry  (0,  a,  0)  and 
(0,  a0,  0)  were  fixed,  the  anticlastic  curvature  of  the  beam  at  that  end  would  give 

w 

rise  to  the  torsion  -^-r^l  («  +  «0)  at  the  free  end.     If  the  fixed  end  be  clamped  in* 


another  way,  the  visible  torsion  due  to  the  anticlastic  curvature  would  be  less  and, 
it  is  clear,  might  even  be  zero.  This  would  happen,  for  example,  if  the  clamping 
were  such  that  it  fixed  that  element  of  the  axis  of  symmetry,  x  =  0,  which,  after 
the  curving  due  to  the  flexure,  remained  in  a  direction  at  right  angles  to  the  plane 
of  flexure.  In  the  practical  case  of  the  beam  which  approximates  both  in  length* 
and  cross-section  to  the  shape  of  the  section  of  a  propeller-blade,  it  is  found  that  the 
anticlastic  curvature  gives  rise  to  a  rotation  about  half  that  due  to  the  shear  torsion 
but  opposite  in  sign. 

It  must  be  borne  in  mind  that  the  anticlastic  torsion  depends  on  the  nature  of 
the  '  fixing  '  and  the  methods  of  clamping  used  in  practice  differ  very  widely  from 
the  theoretical  means  of  fixing  provided  by  the  constants  a,  ft,  y,  a',  ft',  y,  a",  ft"  of 
our  displacement  equations.  These  constants  occur  only  in  the  linear  parts  of  the 
equations  so  that  the  most  we  can  do  with  them  is  to  fix  certain  points  and  certain 

coefficients  of  the  scheme  „  .  '    '  —  {  at  points  on  the  plane  z  =  0.     Remembering  that 

3  (x,  y,  z) 

we  have  found  that  a  =  a",  and  that  there  is  an  equation  for  (ft  —  ft"),  we  see  how  easy 

it  is  to  exhaust  the  stock  of  arbitrary  constants.     Thus,  if  we  try  to  fix  an  element 

.        ,  .  £     ,     ,  .  .  ,  dw      dw      du  I    dv\ 

at  the  point  (x,  y,  0  ),  we  nnd  that  making  u  =  v  =  w  =  0  and  ;—  =  -—  =  —  =0 

dx      cty      dz  \    dx/ 

gives  all  the  arbitrary  constants  definite  values.  Similar  difficulties  arise  with  the 
method  of  loading,  where  the  assumptions  of  the  theory  are  even  more  difficult  to 
approximate  to  in  an  actual  experiment. 

§  10.      The  torsion  of  the  non-curtate  sector  for  various  values  of  the  angle  y. 
1°.     y  very  small. 

To  get  an  approximate  value  of  r  when  y  is  small  we  can  expand  the  trigono- 
metrical functions  and  the  three  series  of  (41)  in  powers  of  y. 

If  &-S  ,-  O=l,2,3], 

;=o  (m  +  2)  my  (m2  —  n2) 


then  SH=  2 

»=o 


m 

ml         m 


-Tit*   T  W"-  «      T      Wl 

•     Z,     -i  nW;,p  -t  jp  I  —  I  «**ap  +  i  -*  2j<  +  i  I  ~  I      f  > 


where  ,,w.!8)+1  =  2  nu9,  nu9+s  =  tf  nu9  +  22""2,  and  nut  =  1.    Further  Tq  =  2  ,    .       .  ,  and  its 

0   \Ali  T  1 ) 

values  are  known  and  tabled,  while  we  ourselves  have  tabled  nuq  from  q  =  4  to  31. 
The  series,  however,  does  not  converge  rapidly  and  can  only  be  used  for  y  small. 

*  For  the  torsion  is  proportional  to  the  length. 
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In  the  latter  case,  omitting  much  troublesome  algebi^a,  we  find  : 

Wa 

T  =  El  kk  +  ^  +  ('°84'033  ~  '^S.OGGyy)  y  -  ('254,632  +  -032,403^)  f 

+  (-027,637-  -198,88277)  y3-  (-084,951-  '130,7767?)  y4  +  etc.].  ..(43  a) 

T/T/yy 

+  '042,017y  -  '262,733y2  -  '022,083y3  -  '052,258y4  -  etc.]  .  .  .(43  6), 


when  we  substitute  the  uniconstant  value  £  for  Poisson's  ratio  77. 
The  limiting  value  as  y-*0  is  therefore 


T=      "x  -28333. 


Wa 

2°.     y  =  3°.     The  approximate  formula  gives  T  =  VJ  x  '28481,  the  accurate  value 

has  the  coefficient  '28504. 
3°.     y  =  6°.     Using  t 

The  accurate  value  of  the  coefficient  is  '28577. 

Wa 

4°.     y  =  12°.    The  approximate  formula  (436)  gives  T  =  ^j  x  '28031,  the  accurate 


3°.     y  =  6°.     Using  the  approximate  formula  (43.6)  we  get  T  —  ^TTJ  x  '28482. 


value  being  "28461.  Clearly,  except  for  very  rough  work,  the  approximate  formula 
cannot  be  used  beyond  about  5°  or  6°. 

In  a  considerable  number  of  cases  the  series  were  summed  and  the  results  are  stated 
below.  This  number  could  easily  be  extended,  but  we  found  the  labour  greater  than 
calculating  by  modern  processes  the  series  to  even  80  or  more  terms.  Great  care  has 
to  be  taken  that  the  right  limits  are  found  for  series  and  terms  which  give  infinite 
values.  These  infinities  always  cancel,  but  the  difference  of  the  two  infinities  may 
provide  a  finite  term.  Interesting  summations  bearing  on  the  theory  of  numbers 
constantly  occur. 

5°.     y  =  45°.     This  case  is  analytically  of  some  interest  for 

»  64  8 

-       ,v=-(i  log,  2  -ITT  +-15886,74780), 
2-  V1 


the  ten-figure  decimal  being  provided  by  Glaisher's  sum  of    2   ,—  =  —  -^  :  see  Proc. 

{  =  0    (4:1  +  0) 

Lond.  Math.  Soc.,  Vol.  vm,  p.  203. 

6°.     y—  90°  —  a  beam  with  right-angled  sector  as  cross-section. 

In  this  case  m  =  2(2i  +  l),  and  the  series  are  all  expressible   in  finite  terms. 
Thus,  omitting  all  the  laborious  algebra, 

64  1 

-  -  4)  -  tan  y  +  y  =  *ff  -  *  8 


»=o 

Substituting  these  in  (41),  we  obtain 
Wa 


*  -  9) 

-T\  =  7     -(16  -Sir  -22  log,  2). 
2)my(m2-!)      180irv 


El'  107r-120loge2 
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7°.     y  =  i80° — a  semi-circular  solid  cylinder  with  flexing  force  parallel  to  the 
diameter  of  the  cross-section.     By  similar  summations  we  find  : 


Wa 


TT 


T  = 


El          7T-8         15 

8°.     y  =  270° — a  solid  circular  cylinder  with  a  quadrant  cut  out  of  it,  the  flexing 
force  being  perpendicular  to  the  bisector  of  this  quadrant : 

Wa      v/2       (4556  -  2352  loge  2  -  5467r)  +  (3068  -  336  log,  2  -  3787r)  -r} 
=  ~EI  X  3T5  ~  2w*-7-8log.2 

9°.     y  =  360° — a  complete  solid  circular  cylinder  with  u  plane  radial  cut  up  to  its 
central  axis,  the  flexing  force  being  perpendicular  to  this  plane  cut. 


We  find  that 


Wa       64  (8 

T  —  -FTF    X 


El      25(97T3-64)' 
Collecting  all  our  results  we  have  the  following  table  : 

TABLE  I.     Torsion  in  a  bent  beam  which  has  for  cross-section  a  non-curtate  sector. 


AngleB  of 
sector 

/Wa, 

Value  of  T  1  j^j  for  any  TJ 

f=i 

*»i 

1=4 

IWa 
r/El 

Total  torsion 

IWa 

Total  torsion 

r  IW<1 

Total  torsion 

Wa/El 

T/EI 

WajEI 

T/EI 

WafEI 

0° 

•066,667+    -866,6671? 

•283,333 

-  -158,333 

•355,556 

-  -188,889 

•500,000 

-  -250,000 

3° 

•070,565+    -857,9117? 

•285,043 

-  -160,043 

•356,535 

-  -189,868 

•499,520 

-  -249,520 

4°-5476 

Interpolated  max.  17  =  ^ 

— 

— 

•356,652 

-  -189,985 

— 

— 

6° 

•073,458+   -849,2577? 

•285,772 

-•160,772 

•356,544 

-  -189,877 

•498,086 

-  -248,086 

6°-9426 

Interpolated  max.  rj  =  ^ 

•285,816 

-•160,816 

— 

— 

— 

— 

12° 

•076,642+    -831,8607? 

•284,607 

-•159,607 

•353,929 

-  -187,262 

•492,572 

-  -242,572 

18° 

•076,772+    -814,0987? 

•280,297 

-•155,297 

•348,138 

-•181,471 

•483,821 

-  -233,821 

24° 

•074,327+    -795,822»? 

•273,283 

-  -145,283 

•339,601 

-•172,934 

•472,238 

-  -222,238 

30° 

•069,762+    -776,9907? 

•264,010 

-•139,010 

•328,759 

-•162,092 

•458,257 

-  -208,257 

45° 

•051,839+    -728,1267? 

•233,870 

-  -108,870 

•294,548 

-•127,881 

•415,902 

-  -165,902 

54° 

•038,645+    -698,4887? 

•213,267 

-  -088,267 

•271,475 

-  -104,808 

•387,889 

-•137,889 

60° 

•029,643+    -679,077r? 

•199,412 

-•074,412 

•256,002 

-  -089,335 

•369,181 

-•119,181 

75° 

•009,098+    -633,7077? 

•167,524 

-  -042,524 

•220,333 

-  -053,666 

•325,951 

-  -075,951 

90° 

-•004,459+    -596,1547? 

•144,580 

-  -019,580 

•194,259 

-  -027,592 

•293,618 

-  -043,618 

104°-8863 

Interpolated  min.  T?  =  j 

•136,646 

-•011,646 

— 

— 

— 

— 

106°-1906 

)!                          )'                *?=J 

— 

— 

•184,139 

-  -017,472 

— 

— 

.108°-3252 

..              »         i?  =  J 

— 

— 

— 

— 

•278,781 

-  -028,781 

120° 

•006,222+    -558,6557? 

•145,886 

-  -020,886 

•192,440 

-  -025,773 

•285,549 

-  -035,549 

150° 

•082,335+    -584,3807? 

•228,430 

-  -103,430 

•277,128 

-  -100,461 

•374,525 

-  -124,525 

180° 

•224,047+    -672,1417? 

•392,082 

-  -267,082 

•448,094 

-•281,427 

•560,117 

-•310,117 

210° 

•407,389+    -798,7377? 

•607,073 

-  -482,073 

•673,635 

-  -506,968 

•806,758 

-  -556,758 

240° 

•597,136+    -929,6507? 

•829,549 

-  -704,549 

•907,019 

-  -740,352 

1-061,961 

-•811,961 

270° 

•755,377+l-028,203r? 

1-012,428  j  -  -887,428    1-098,111 

-  -831,444 

1-269,478 

-1-019,478 

300° 

•852,715  +  1-066,4257? 

1-119,321 

-•994,321 

1-208,190 

-1-041,423 

1-385,928 

-1-135,928 

319°-4474 

Interpolated  max.  T?  =  | 

1-138,590 

-1-013,590 

— 

— 

— 

— 

318°-2777 

»              »         ??  =  i 

— 

— 

1-226,331 

-1-059,664 

— 

— 

316°-3202 

,.          i)      ?=i 

— 

—                —  . 

— 

1-402,199 

-1-152,199 

330° 

•875,056+1-031,7317? 

1-132,989 

-1-007,989  :  1-218,967 

-1-052,300 

1-390,922 

-1-140,922 

360° 

•824,927+    -928,0437? 

1-056,938 

-•931,938    1-134,275 

-  -967,608 

1-288,948 

-1-038,948 

The  total  torsion  =  anticlastic  torsion  minus  shearing  torsion,  and  it  tends  to  move  the  angle  edge  of 
prism  in  the  opposite  direction  to  the  flexing  force.  Since  the  total  torsion  is  always  negative  the  angle 
edge  of  prism  will  always  move  in  direction  of  flexing  force. 
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The  curve  in  Fig.  3  shows  how  the  nature  of  the  torsion  changes  with  y.  It 
provides  a  kind  of  physical  measure  of  the  asymmetry  of  the  circular  sector,  for, 
as  we  have  seen,  the  effect  is  wholly  due  to  the  distribution  of  shearing  forces  over 
an  asymmetrical  section.  Remembering  the  direction  of  the  torsion  as  given  by 
equation  (42)  and  Fig.  2,  we  see  that  in  the  case  of  a  wedge-shaped  sector  of  small 
angle,  such  as  the  blade  of  a  sword  held  horizontally*,  the  torsion  will  be  such  that 
the  back  twists  up  with  a  downward  flexing  force.  A  first  maximum  is  reached 
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Angle  of  sector  (y) 

FIG.  3.     Torsion  clue  to  Shearing  Forces  in  a  flexed  beam  the  cross-section  of 
which  is  a  non-curtate  circular  sector. 

between  4°  and  7°  according  to  the  value  of  Poisson's  ratio.  As  we  increase  the 
angle  of  the  sector  the  torsion  diminishes  in  magnitude  until  with  y  equal  to  about 
105°  the  value  of  T  is  a  minimum.  From  consideration  of  asymmetry  this  is  what  we 
might  expect,  for  with  y  =  105°  we  have  from  various  points  of  view  the  sector  which 
most  nearly  approaches  symmetry.  With  angles  larger  than  this  angle  of  approximate 
symmetry  the  torsion  of  the  beam  increases  again.  With  the  increasing  departure 

*  A  rough  experiment  was  made  with  a  carving-knife  to  which  a  long  light  pointer  was  attached 
to  show  the  torsion  when  a  weight  was  hung  on  the  end  of  the  blade.  The  torsion  was  of  course 
extremely  small,  but  was  seen  to  be  in  the  direction  given  by  the  theory. 
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from  symmetry  as  y  becomes  larger,  the  amount  of  torsion  continues  to  increase  until 
a  second  maximum  is  reached  about  316°  to  319°,  from  which  value  it  falls  slightly 
until  y  —  2ir  which  is  the  physical  limit  to  the  application  of  the  formula. 

§11.     Numerical  examples  of  the  non-curtate  sector. 
The  total  torsion  in  a  beam  of  length  I  is  given  by 


rl  = 


Wai 
El 


x  the  constant  tabled  in  Table  I 


Wl*      3a     A 

~  o  ZT rr  x  ~JT  x  the  tabled  constant 
o-tLl        I 

=  S  x  -TT  x  the  tabled  constant, 

if  8  be  the  deflection  of  the  centroid  of  the  end  section.  This  is  on  the  assumption 
that  the  length  I  is  so  great  that  we  can  neglect  the  linear  term  sj  (see  equation  (50) 
infra)  which  comes  into  the  formula  for  the  deflection  owing  to  the  shearing  forces, 
an  assumption  which  we  shall  investigate  later  in  the  paper*. 

TABLE  II.     The  deflection  and  torsion  of  steel  beams  of  circular  sector 

cross-section  with  a  load  W  Ibs. ,  for  varying  size  of  angle. 
Length  of  beam  =  20  inches,  breadth  =  1  inch,  stretch-modulus  =  30,000,000  Ibs.  per  square  inch. 


Angle  of 
sector  =  7 

WP 
Deflection  8  =         inches 
o£i 

Shear  torsion  =rl  degrees 

Shear  torsion  per 
deflection  in  degrees 
per  inch 

r 

IF  x  802  -53 

fFx97-9356 

0°-122 

3° 

x    29-73 

x    3-6416 

0°-122 

6" 

x      3-717 

x    0-45646 

0°-123 

12° 

x      0-465 

x    0-05686 

0°-122 

18° 

x      0-138 

x    0-01662 

0°-120 

24° 

x      0-0586 

x    0-00688 

0°-117 

30° 

x      0-0301 

x    0-003,415 

0°-113 

45° 

x      0-00908 

x    0-000,9123 

0°-100 

60° 

x      0-00393 

x    0-000,3368 

0°-086 

90° 

x      0-001,246 

x    0-0000,7741 

0°-062 

120° 

x      0-000,5789 

x    0-0000,3629 

0°-063 

150° 

x      0-000,3357 

x    0-0000,3295 

0°-098 

180° 

x      0-000,2264 

x    0-0000,3815 

0°-168 

210° 

x      0-000,1707 

x    0-0000,4453 

0°-261 

240° 

x      0-000,1407 

x    0-0000,5016 

0°-356 

270° 

x      0-000,1245 

x    0-0000,5417 

0°-435 

300° 

x      0-000,1165 

x    0-0000,5604 

0°-481 

330° 

x      0-000,1136 

x    0-0000,5531 

0°-487 

360° 

x      0-000,1132 

x    0-0000,5141 

0°-454 

NOTE.    A  positive  torsion  signifies  a  torsion  in  the  opposite  direction  to  the  arrow  in  Fig.  2,  p.  20. 

In  Table  II  we  have  as  a  particular  case  the  data  for  a  steel  beam  of  length 

20  inches  and  radius  1  inch,  the  load  applied  being  W.    Taking  E  to  be  30,000,000  Ibs. 

a4 
per  square  inch  and  calculating  /  from  the  formula  —  (y  —  sin  y)  we  get  the  deflection 

o 

*  See  §  14. 
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and   the  torsion  of  the  beam.     For  another  beam  of  length  I'   inches,   of  radius 

a'  inches,  and  stretch-modulus  E',  the  deflection  is  got  by  multiplying  the  deflec- 

l'3 
tions  of  Table  II  by  3,750  ,-,,  ,4  and  the  total  torsion  by  multiplying  the  tabled  values 

by  1,500,000  --. 


§  12.     Distribution  of  vertical  shear  along  the  central  radius  of  a  cross-section. 

To  get  some  idea  of  the  distribution  of  shears  which  gives  rise  to  the  torsion,  we 
will  consider  the  distribution  of  the  vertical  shear  along  the  central  radius  of  any 
cross-section,  that  is  the  shear  zx  along  the  line  x  =  Q  in  any  cross-section*.  We  shal] 
confine  ourselves  to  the  case  of  the  complete  sector. 

From  equations  (16)  and  (40), 


W 

TV-^ 

• 


~  m 

dx\_  2  cos  y 


*  sin  20 


Along  x  —  0,  or  6  —  0  we  have 


\cosy 


8r(-l)'     1 
'2-m-2 


El  y      \2(m2-F)      2(m2-9)/  a 

when  we  put  rj  =  ^  and  insert  the  expression  for  Am, 


El  \\'    El    a  \cosy 

v  »     I  \  / 

Wa 


.j/ 

t  =  0\ 


i    •         if 
where  m  —  (2i  +  1 )  - . 

7 

*  Shears  other  than  in  the  vertical  direction  will  of  course  have  appreciable  and  in  certain  cases  may 
be  of  great  effect  on  the  value  of  T,  but  the  labour  of  calculation  would  be  considerable  except  in  the  case 
we  have  taken,  and  that,  after  all,  is  the  one  likely  to  be  most  instructive.  An  important  result  also  flows 
from  this  central  radius  shear :  it  enables  us  to  determine  (see  p.  44)  the  shear  at  the  centroid  which 
provides  the  linear  term  of  the  deflection  with  Saint-Venant's  form  of  the  terminal  fixing. 
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The  form  of  equation  (45)  allows  us  immediately  to  compare  the  distribution 
of  shear  in  beams  of  various  angles,  but  of  constant  length  and  radius,  when  the 

1  Wl3 
deflection  ~  VT   is  kept  constant.     This  is  the  suitable  method  of  comparison,  and 

O    fj  1 

I       TT7     °\ 

in  what  follows  we  will  only  consider  zx   (  fi   UT  )  ,  that  is  the  expression  in  the  double 


.,  ,  , 

brackets,  it  bem£  proportional  to 


the  vertical  shear  on  central  axis 
—  :  —  -  . 

deflection 


It  is  useless  to  attempt  to  see  the  nature  of  the  distribution  of  shear  from  the 
equation  itself,  and  we  have  immediately  to  resort  to  numerical  calculation. 

In  Table  III  below  are  given  the  values  of  zx/ifi  ^,T  )  for  every  tenth  of  the 

\         -*-/-*.     / 

radius  of  the  sector  and  for  sectors  whose  angles  are  5°,  12°,  30°,  45°,  60°,  72°,  90°, 
120°,  150°,  180°,  270°,  360°. 

For  all  these  values  of  y  the  series  involved  in  equation  (45)  rnay  each  be  summed 
in  finite  terms.     For  let  us  consider  the  series 


j»om--jr\ay 

a  series  which  for  p=  1,  2,  3  gives  the  three  series  concerned. 
1 


m2  —  p2  \a 


1     f/i/X11"1          fa  /y\m~p~l 

(i)      2        ^  (- 

2p\_\a/          Jo\a/ 


~p~1 


y 

a 

oa 


m+p-i 


a 


=  JLT/y 

2p[\a 


'?• 


jzy -nil)-:      '(-t,*tifi\. 

W       JOi-oW  \a/J 

the  interchange  of  the  integral  and  summation  signs  being  clearly  justified, 


2p 


\a 


1+1* 

\a 


1+  * 

\a 


(4G), 


and  both  of  these  integrands  being  rational  functions  of     /- ,  ivhen        is  an  integer 

\'  Cl  y 

(as  in  all  the  cases  considered),  we  must  be  able  to  express  the  series  in  finite  terms 
by  evaluating  these  integrals.     In  general  they  are  incomplete  B-functions. 
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For  the  sectors  whose  angles  are  5°,  12°,  30°,  45°,  60°,  72°  and  120°,  it  is  easier  to 
use  the  series  as  they  stand  than  to  find  the  sums  in  finite  terms  and  calculate  from 
them.  For  these  cases  the  equation  (45)  becomes  : 


^  //    W<n~\  v  fij\-  /?/\35  ///\107 

for     5°,      zx/[u.^=rf)=  +  '001091^  +  '000959  K     -'002668K      +  '000283K       +..., 
/  V   El  /  a  \al  \a/  \a] 

•-^  //    W<y*\  ij  /?/\2  /ij\u 

for  12°,     zx     it-p"-   =  +  '006358^-  +  '005713  (y-\  -  '007971  (^ 
/  V    LI  I  a  \a/  \a/ 


/?/\44  iv\'4  ii/\m 

+  '000701  K)    -'000247   ^     +  '000128  K      +..., 
\aj  \aj  \a/ 

for  30°,     zx/(u.~~]=  +'040842^  +  '045753  (^  -  '044065  (%-} 
/  Y    El  I  a  \a/  \aj 

+  '001400  C^)    -'000426  AM    -..., 
\a/  \a/ 

^^  //    Wa2\  ?/  /?v\2  /w\3 

for  45°,     zx   (u.-trr  }=  +'096872 ^  +  '176777    ^     -'176109   ^ 
/  \    El  /  a  \a]  \a] 

+  •001295^)    -'000210^)   +'000073^)+..., 
\a/  \a/  \a/ 

-^  //    Wo?\  ij  /V-2  /7y\2         11 

for  60°,      zaJ/flt-^r  1-  +  '199412^  -  '033411    ^  )  -  '4760556    "    log]0^  +  "000832 
/  \    El  /  a  \a]  \al      '  °a 

+  -000286  (^)    -'000221^)   +'000152^)    - '000108  f'^   +.. 
\a/  \a/  \a/  \a/ 


for  72°,      9/u  —   -  +  -38699  y-  -  '452254  ($ 

El  /  a  a 


ly  /  y  /v\e  /y\n  /W\16  \ 

+  \y-  '289781  ^  +  '000780  (*-  +'000706^  -  '000481  "  +...  , 
\/  a\      a        \a/         \a/         \a)     / 

•^  1  1     ~WrP\  11  /i/\2 

for  120°,    zx   (p  j~}=-  '437658  ^--'187500  ^1 
/  \^  El  I  a  \aj 

fa  I  /w\3  /v\e  /w\9 

+     fl    1-071134  +  -022715   £)  -'003120  ^    +'001121   ^ 

v*.\  W  \fll  W 

/i/\12  /?/\15  /W\18  /V/V1 

-  -000579    ^     +'000355    '^     -'000241  (^     +  '000209  (J  ) 

\aj  \a]  \aj  \aj 


•000134  f^V+  -000105  f^V-  -000084  f^V  + '000070  f'- 
\a/  \a/  \a/  \a 


-•000059    ^     +'000050    ^      - '000043   ^ 
\aj  \a' 


For  the  sectors  whose  angles  are  90°,  180°,  360°  the  series  are  too  slowly 
convergent  for  calculation,  but  the  sums  of  the  series  are  without  much  difficulty 
obtainable  by  the  method  used  to  obtain  equation  (46). 
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For  y  =  90°: 

Wa*\ 

zx  i 


29 


l(y. 

4\a 


4      4 


1  /y  V 


f272tan"      (y 


4  \a/         4  \aj  /     4  72 


1  /« 

+  7  (- 
4  \a     va 


1-f* 

a 


!> 

1  .        tt- 

— i-tan"1 - 

272  /?/ 

u 


.(47). 


Instead  of  this  complicated  expression  the  series  formula  of  equation  (45)  may  be 

I/ 

used  directly  for  low  values  of  - .     Actually,  in  the  calculations  for  the  table  given 
on  p.  31,  this  method  was  used  as  a  check  on  the  accuracy  of  the  work. 

For  7  =  180°,  the  case  of  a  semi-cylinder  flexed  in  the  plane  parallel  to  its  plane 
side : 


zx 


8      8 


a 


a 


a 


(48). 


For      =  270" 


zx 


a 


,48-50  K 
a 


1  /V 

o     - 

3  \a 


where 
and 
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For  y  =  360°,  the  case  of  a  complete  cylinder  with  a  radial  cut  perpendicular  to 
its  plane  of  flexure  : 

Wcf\_  _     I  jWa\\ 


•-F  ly-(-1 

47T  |_V  a  \      3 


y    >r\±    -i   A/1 

-  —  -       tan  '    /- 

a     a    J  \/  «J 


3      a       5  \a/          5  \aj         3\al         \a/ 

?1 (50). 


\a/         \aj 

In  the  cases  of  y=60°,  y  =  90°,  y=  1803,  m  takes  the  values  3,  2,  1  respectively 
and  consequently  there  arise  infinite  terms  in  the  series  of  equation  (45).  There 
are  however  corresponding  trigonometrical  infinities  of  the  same  order  and  the 
limits  of  the  whole  expression  are,  as  they  must  be,  finite.  For  y>  180°  the  shear 

becomes  infinite  at  -  =  0.     This  is  the  case  of  the  re-entrant  angle  referred  to  later 
a 

on  in  this  section. 

The  accompanying  Table  III  gives  the  collected  numerical  results  and  in  Fig.  4 
are  diagrams  each  representing  the  distribution  of  shear  along  the  central  radius 
for  different  sectorial  angles.  It  would  be  impossible  to  show  adequately  in  one 
diagram  the  comparative  distributions  of  stress  for  the  different  sizes  of  the  sectorial 
angle.  It  is  to  be  remembered  in  considering  the  table  and  the  diagrams  that,  in 
accordance  with  the  axes  of  reference  which  we  have  taken,  a  shear  of  the  material 
on  the  fixed  side  on  material  on  the  loaded  side  of  the  cross- section  is  positive  when 
it  acts  opposite  to  the  direction  of  load. 

Our  diagrams  provide  the  excess  of  actual  shear  over  mean  shear  in  terms  of 
mean  shear.  They  indicate  the  point  at  which  the  shear  is  exactly  equal  to  the  mean 
shear  and  the  value  of  the  shear  at  the  centroid.  When  y  =  0  we  have,  if  A  stand  for 
area  of  section  and  WjA  for  mean  shear, 

zz/(  W/A)  =  1-36^+  1-2  (2 
a  \a 

or  the  limiting  case  is  a  parabola.  The  ratio  runs  from  0  to  2  "5 6,  or  the  excess  on 
mean  shear  from  —  1  to  1'56.  The  curve  is  given  in  Fig.  4  (i)  of  this  series.  Almost 
at  the  half  radius  ('508)  the  shear  is  equal  to  the  mean  shear,  but  on  the  right  of  this 
point,  i.e.  towards  the  thick  end,  the  excesses  on  the  mean  shear  are  considerably 
greater  than  the  defects  on  the  left,  and  we  get  light  on  why  the  thick  end  rotates 
upwards.  We  may  look  upon  this  curve  as  the  parabola  plus  the  vertical  line  through 
the  point  1  '0.  This  explains  how  the  apex  curve  comes  away  as  in  Fig.  4  (ii)  of  the 
series,  when  y  =  5°,  and  the  vertical  line  bends  round  so  that  we  actually  get  negative 
central  shears  at  the  butt  end  in  sector  prisms  of  very  small  angle.  A  horizontal  line 
through  —  TO  of  vertical  scale  indicates  in  all  these  diagrams  the  extent  of  the 
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FIG.  4.    Diagrams  showing  the  Distribution  of  Vertical  Shear  along  the  central  radius  of 
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positive  and  negative  central  radius  shears.  Passing  to  Fig.  4  (iii)  we  find  for  y  =  12° 
the  defect  on  the  mean  shear  still  negative  but  very  slight  ( —  '023),  the  position  of 
maximum  shear  is  passing  away  from  the  butt  end,  and  the  position  of  mean  shear  is 
very  slightly  to  the  left  of  the  mid-point  of  the  radius.  These  features  are  still  more 
marked  in  Fig.  4  (iv),  y  —  30°,  where  we  have  a  considerable  excess  over  the  mean  shear 
at  the  butt  end,  the  position  of  mean  shear  is  visibly  separated  from  the  mid-point 
and  the  shear  at  the  butt  end  still  increasing.  The  process  is  continued  still  further  in 
Fig.  4  (v),  y  =  45°,  but  the  butt  end  shear  in  terms  of  mean  shear  has  now  reached  a 
maximum  and  with  Fig.  4  (vi),  y  =  60°,  has  begun  to  fall.  The  shear  at  the  centroid  is 
now  also  falling ;  and  the  point  of  mean  shear  is  moving  rapidly  towards  the  apex. 
With  Fig.  4  (vii),  y  =  72°,  the  tendency  of  the  curve  to  flatten  out  horizontally  is  con- 
spicuous and  this  is  accompanied  by  a  fall  in  the  butt  end  shear  and  the  centroid  shear 
with  a  still  further  motion  of  the  point  of  mean  shear  towards  the  apex.  In  Fig.  4  (viii) 
and  (ix)  for  y  =  90°  and  y=120°  respectively  these  features  are  emphasised.  In 
Fig.  4  (x)  for  y=  150°  we  are  approaching  a  climax.  The  butt  end  shear  is  now  less 
than  the  mean  shear,  the  position  of  mean  almost  up  to  the  apex  and  the  flattening 
of  the  top  of  the  curve  now  stretches  across  from  0'2  to  0'5  almost  of  the  radius. 
A  section  between  y  =  150°  and  y  =  180°  would  show  the  curve  with  a  pointed  apex 
like  Fig.  4  (ii)  reversed  right  and  left,  and  this  enables  us  to  interpret  Fig.  4  (xi)  which 
is  almost  a  straight  line.  For  y  =  180°  the  first  point  of  mean  shear  has  reached  the 
apex,  the  zero  shear  has  ceased  and  we  have  a  finite  excess  over  the  mean  shear  at 
the  apex ;  the  curve  indeed  should  be  considered  as  built-up  of  the  nearly  straight 
sloping  line  and  the  vertical  through  O'O.  The  second  point  of  mean  shear  has  reached 
the  neighbourhood  of  0'8  and  we  have  increasing  defects  on  mean  shear  at  the  butt 
end.  Fig.  4  (xii),  y  =  270°,  shows  the  revolution  which  follows  on  the  re-entrant  angle. 
The  shear  at  the  apex  is  now  infinite  ;  the  shear  at  the  butt  end  is  now  negative  and 
probably  approaches  a  maximum  somewhere  about  300°.  The  position  of  mean  shear 
has  also  reached  somewhat  about  its  maximum  in  motion  towards  the  right.  In 
Fig:  4  (xiii),  y  =  360°,  or  the  solid  cylinder  with  a  cut  to  the  centre,  the  butt  end 
shear  is  again  positive  and  less  than  the  mean  shear,  the  position  of  mean  shear  has 
returned  towards  the  butt  end,  and  the  big  shears  are  all  towards  the  left,  the  exact 
opposite  to  what  we  noted  in  the  sectors  of  small  angle  first  discussed,  but  at  the- 
same  time  they  are  all  to  the  right  of  the  centroid,  so  that  their  total  turning  moment 
is  big  and  the  result  a  large  couple  turning  the  butt  end  upwards. 

The  minimum  torsional  effect  is  shown  in  Table  I  to  occur  in  a  sector  of  about 
105°.  This  is  roughly  indicated  in  Fig.  4  (viii)  and  (ix)  where  the  couple  of  the 
central  radius  shears  at  least,  as  measured  by  the  moments  of  the  areas  above  the 
horizontal  line  through  —  1 -0  about  the  centroid,  appear  mostly  closely  balanced. 

Of  course  the  reader  must  bear  in  mind  that  these  central  radius  shears  are  not 
the  whole  shearing  system,  but  the  central  radius  is  the  locus  probably  of  maximum 
shear  as  in  the  case  of  the  flexure  of  symmetrical  systems. 
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§  13.      The  torsion  of  the  curtate  sector  for  particular  cases. 

The  formulae  for  the  curtate  sector  are  very  much  more  complicated  than  the 
formulae  for  the  complete  sector  and  require  a  large  amount  of  additional  arithmetic. 
This  is  however  much  lessened  if  the  case  of  the  complete  sector  for  the  angle  in 
question  has  been  worked  out,  for  the  calculation  then  consists  of  multiplying  all  the 
quantities  concerned  by  appropriate  factors.  These  factors,  also,  converge  to  unity 
with  great  rapidity  if  the  ratio  of  the  inner  radius  a0  to  the  outer  radius  a  is  not  too 
large  and  if  the  value  of  y  is  fairly  small  so  that  m  is  large.  This  is  the  case,  for 
instance,  when  y  =  12°,  which  is  given  below. 

For  the  limiting  case  of  «„-»«,  we  can  expand  all  the  factors  of  the  form  ap  —  aap 


___  ft 


___ 

in  powers  of  e  =  -      -  ,  and  we  obtain  for  limiting  value 

Wn  Wn 

(84'5G  +  86707?)  =  -     r  x  106-24,  if  T?  = 


1°.     A  thin  wedge-shaped  curtate  sector. 

The  angle  chosen  was  12°,  and  the  results  of  the  calculations  are  given  below  in 

/  Wa 
Table  IV  and  graphically  in  Figs.  5  and  6.     The  column  headed  i/»,  =  T/-pT  gives  the 

constants  as  directly  calculated,  and  the  following  column  headed  &  =  T7  —  v,       -  is 


deduced  from  the  previous  one.    The  first  expression  gives  the  ratio  of  total  torsion  (rl) 

/8      1  WF       h  3 


„     ,.  .   . 

to  deflection  per  unit  length  I  -,-  =  -  jsj\  =  j  x  —  -  -  —  TTLVii  wnere  »  1S  the  maximum 

height  of  section  and  i/^  the  number  in  the  second  column  of  Table  IV.    It  shows  that 
for  constant  height  of  section  and  length  of  beam  the  torsion  at  first  decreases  slowly 

TABLE  IV.      The  torsion  per  unit  length  in  a  sector  of  central  angle  12° 
and  external  radius  a  for  varying  degrees  of  curtateness. 


Torsion  per  unit  length 

Internal  radius      o0 

External  radius      a 

IWa 

fc"r/iz 

*-Pa* 

o-o 

•2846 

•2846 

•3 

•2486 

•3729 

•6 

•220.r> 

•6614 

•75 

•2219 

•8875 

•8 

•2248 

1-1238 

•9 

•2386 

2-3861 

•95 

•26S6 

5-3722 

•98 

•4528 

22-6390 

1-00 

106-24 

oo 

NOTE.     A  positive  torsion  signifies  a  torsion  in  the  opposite  direction  to  the  arrow  in  Fiji;.  2,  p.  20. 
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and  then,  after  a0  =  -9a,  increases  rapidly  with  the  squatness.  The  second  expression 
indicates  that  the  same  ratio  may  be  represented  by  — ~ ?  x  3i|/.,,  where  i/f.,  is  the 

number  in  the  third  column  of  Table  IV,  and  the  values  of  \jjz  show  that  for  constant 
breadth  of  section  and  length  of  beam  the  ratio  of  total  torsion  to  deflection  per  unit 
length  increases  at  first  slowly  and  then  very  rapidly  with  the  squatness. 


4- 


Torsion  measured  on 
Breadth  of  Section 


Torsion  measured  on 
Height  of  Section 


FIG.  5. 


FIG.  6. 


00 


O'O 


FIG.  5.    The  torsion  of  a  flexed  beam  which  has  as  cross-section  a  thin  (angle  =  1 2°)  wedge-shaped 
sector  with  varying  curtateness,  the  torsion  being  measured  on  the  external  radius. 

FIG.  6.     Torsion  of  the  same  beam  as  in  Fig.  5,  the  torsion  being  measured  on  the 
central  breadth,  not  the  external  radius. 
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Numerical  example. 

As  in  the  case  of  the  non-curtate  sector,  we  shall  now  consider  the  special  example 
of  a  beam  of  steel  20  inches  long  having  as  cross-section  a  curtate  sector  of  angle  12° 
and  outer  radius  1  inch  (the  outer  arc  is  therefore  in  all  cases  '209...  inches).  The 
values  of  the  deflection  and  torsion  when  a  load  W  is  applied  are  given  in  Table  V 
and  the  torsion  per  deflection  for  varying  curtateness  is  shown  in  Fig.  6,  the  scale  in 


degrees  per  inch  being  obtained  from  the  scale  of 
by  the  factor  "430. 


Wa 


there  given  by  multiplying 


TABLE  V.     The  deflection  and  torsion  of  a  steel  beam  of 
curtate  circular  sector  cross-section. 

Length  of  beam  =  20  inches. 

Outer  arc  =  0-209  inches. 

Angle  of  sector  =  12°. 

Stretch-modulus  =  30,000,000  Ibs.  per  square  inch. 


Internal  radius      a0 

Deflection 
1  Wl3  . 
=  3£lmCheS 

Total  torsion 
=  T  I  degrees 

Torsion  (in  degrees) 
per  inch  of 
deflection 

External  radius      a 

0 

Wx    -465 

Wx  -0569° 

0°-1223 

•3 

x    -471 

x  -0503° 

0°-1068 

•6 

x    -580 

x  -0549° 

0°-0948 

•75 

x    -681 

x  -0649° 

0°-0954 

•8 

x    -788 

x  -0762° 

0°-0966 

•9 

x  1-353 

x-1389° 

0°-1027 

•95 

x  2-509 

x  -2896° 

0°-1154 

1-0 

X  OO 

x  QO 

45°  -68 

NOTE.     A  positive  torsion  signifies  a  torsion  in  the  opposite  direction  to  the  arrow  in  Fig.  2,  p.  20. 

2°.     A  cylindrical  tube  with  a  longitudinal  cut. 

If  we  take  y  =  2ir  we  get  a  sector  which  is  a  full  circular  ring  except  that  it  has 
a  cut  along  one  radius.  The  beam  is  then  really  a  cylindrical  tube  with  a  cut  along 
one  side.  This  case,  especially  if  the  tube  be  thin  (e.g.  as  below  aa/a  =  -9),  is  of 
peculiar  interest  as  it  represents  the  extreme  of  asymmetry  in  a  circular  sector.  The 
method  of  approximation  for  a0-*a  mentioned  at  the  beginning  of  this  section  having 

failed,  we  have  to  use  direct  calculation  and,  both  y  and  the  ratio  —  being  large,  the 

ct 

calculation  of  the  series  is  extremely  heavy.  The  factor  of  T  in  equation  (39)  turns 
out  to  be  the  difference  of  two  nearly  equal  numbers  and  to  get  three  figure  accuracy 
in  the  final  result  it  was  found  necessary  to  calculate  a  large  number  of  terms  of  the 
series  to  10  decimal  places. 

If  we  take  a0  =  '9a  we  get 


Wri 


r= 


520-6 


Wn 


650-75,  for  r,  = 


as  compared  with 


r-45x  1-057 
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for  a  non-curtate  sector  of  the  same  angle.  This  is  an  extraordinarily  large  value  for 
T,  but  it  is  well  known  that  a  hollow  rod  loses  much  of  its  power  to  resist  torsion 
when  it  has  a  cut  down  one  side  and  it  might  therefore  be  anticipated  that  the  couple 
due  to  the  great  asymmetry  of  the  flexural  shears  would  bring  about  a  correspondingly 
large  torsion. 

We  can  illustrate  this  result  by  taking  the  value  of  T  for  e  =  (a  —  «0)/«j  very  small. 
The  deduction  of  this  value  is  of  much  interest  as  it  involves  the  summation  of  a  number 
of  noteworthy  series  which,  however,  cannot  be  considered  here  ;  ultimately  we  find  : 

Wa  47124(1  +•>?)_  Wa  _  Wa  5-8905 
~  ~~~  "~~~~~ 


T~ 


EI 


indicating  the  high  value  of  T  as  the  thickness  diminishes  indefinitely*. 

Experimental  Verification. 

With  such  a  large  value  of  T  it  is  easily  possible  to  see  the  twisting  of  such 
a  beam  under  flexure,  and  some  rough  preliminary  experiments  were  made  with  a 
piece  of  piping  used  as  a  casing  to  electric  wiring.  This  is  manufactured  with 
the  required  slit.  One  end  was  fixed  in  the  chuck  of  a  turning  lathe,  a  solid 
boss  having  been  put  inside  the  tube  at  that  end  to  prevent  buckling.  The  load  was 
applied  by  a  stirrup  arrangement  fastened  to  the  other  end.  The  deflection  was 
measured  by  taking  the  vertical  heights  of  the  bottom  of  tube  from  the  lathe-bed  or 
from  the  slide-rest.  The  torsion  of  the  load-end  relative  to  the  fixed  end  was  measured 
with  the  aid  of  a  light  pointer  fastened  rigidly  at  the  load-end  moving  over  a  plate 
fastened  rigidly  to  the  fixed  end.  The  effective  radius  of  the  pointer  was  10  inches 
and  the  torsion  was  measured  on  a  12"  length  of  the  15"  pipe  so  that  with  the  largest 
observed  torsion  of  '159  radians  or  9°  the  pointer  moved  through  l-59  inches,  and  the 
real  torsion  was  ^f  .  159  =  '199  radians. 

The  dimensions  of  the  tube  were  :  Outer  radius  =  0'50  inches,  inner  radius  =  0'45 
inches,  length  (excluding  the  part  used  in  fixing)  =  15  inches. 

TABLE  VI.     Experiments  with  a  piece  of  split  piping. 
Deflection  and  torsion  with  various  loads. 


Load  in 
lbs.  =  JT 

Experimental 

Theoretical  (£  =  27,000,000) 

Average  deflec- 
tion at  end  in 
inches  =  5 

Average  torsion  =  rl 

Deflection  in 
inches  (linear 
and  cubic  terms) 

Torsion  in  radians 

In  radians 

In  degrees 

«0  =  -45 

n0  =  -425 

0 

•000 

•000 

0°    (I' 

•000 

•000 

•000 

3-7 

•048 

•015              0°  52' 

•037 

•041 

•017 

13-7 

•131 

•059 

3°  23' 

•138 

•151 

•063 

23-7 

•215 

•105 

6°    1' 

•239 

•2G2 

•108 

33-7 

•305 

•152 

8°  43' 

•339 

•372 

•153 

43-7 

•387 

•199 

11°  24' 

•440 

•488 

•199 

*  If  e  be  put  O'l  we  have  r  =  ^y589'05,  which  is  of  the  same  order  as  the  true  value  above,  but 
indicates  that  r\j  is  not  small  enough  for  us  to  neglect  the  next  term  in  e. 
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In  the  first  experiment  the  deflection  and  torsion  were  measured  only  at  the  end 
of  the  beam.     The  results  are  given  in  Table  VI  and  represented  in  Fig.  7. 


o 
H 


T3 
03 

Q 


,  Tors'0" 


0  5  10          15          20          25         30          35          40         45 

Load  in  Ibs.,  W 
FIG.  7.     Deflection  and  Torsion  with  varying  load  of  a  split  pipe*. 


o 


5  10 

Distance  along  beam  in  inches,  z 
FIG.  8.     Curve  of  Flexure  for  split  pipe. 


15 


It  is  evident  that  both  the  deflection  and  the  torsion  are  proportional  to  W,  but 

there  is  a  great  divergence  between  the  observed  values  of  8  and  the  values  that 

1  Wl3 
are  given  by  direct  use  of  the   Eulerian  expression  i^-prf  for  S,   E  being  taken 

to  be  27, 000, 000  f  Ibs.  per  square  inch.  This  leads  us  to  question  whether  the  linear 
term  s0z  in  the  expression  for  the  displacement  is  negligible  as  it  is  often  taken  to 
be  in  beams  of  moderate  length.  As  we  shall  see  the  theory  does  not  give  much 

*  The  deflections  were  taken  on  a  15",  the  torsions  on  a  12"  length,  so  that  the  latter  need  to  be 
multiplied  by   1'25. 

t  This  choice  is  based  on  calculations  given  later  on  in  this  section. 


40 


YOUNG,  ELDERTON  AND  PEARSON 


help  because  the  nature  of  the  clamping  is  so  uncertain,  but  an  experimental  test  is 
easy  and  obvious.  If  we  measure  the  deflection  at  various  points  along  the  beam  we 
ought  to  be  able  to  pick  out  the  linear  term  if  it  be  not  negligibly  small.  Observations 
were  therefore  taken  of  the  deflection  of  the  beam  with  a  load  of  437  Ibs.  at  three 
points  along  the  beam.  The  means  of  several  trials  are  as  follows. 

TABLE  VII.     Experiments  with  a  piece  of  split  piping. 
Observations  to  determine  curve  of  flexure. 


Distance  along  beam  (z)  in  inches: 

0 

4i 

9 

15 

Average  deflection  8  in  inches 

0 

•070 

•166 

•313* 

Approximation  of  equation  (51)  

0 

•068 

•169 

•313 

These  are  plotted  in  Fig.  8  on  p.  39.     If  we  assume  the  theoretical  formula 

W 
S,  =  sf+  jjj  (ife'-iz3) (51), 

and  substitute  the  values  of  W,  I,  1,  we  find  by  the  method  of  least  squares  that 

«„  =  '014, 

and  E  =  27,000,000, 

giving  8Z  =  '01 4z  +  '00072z2  -  "00001 6z3. 

This  value  of  E  is  quite  near  what  might  be  expected  for  this  grade  of  mild  steel. 
The  large  value  of  s0  appears  at  first  extraordinary  but  it  is  even  less  than  some 
theoretical  forms  of  clamping  '  provide ;  with  it  we  should  have  to  take  a  length 
of  over  65  inches  before  the  linear  term  was  less  than  10  °/0  °f  the  cubic  term  in 
the  expression  for  the  end-displacement.  Exactly  the  same  experiment  was  made 
on  the  pipe  rotated  axially  through  90°,  so  that 
instead  of  having  the  cut  to  the  side  as  in  (a)  of 
Fig.  9  it  was  vertically  upwards  as  in  (b),  and 
the  least  square  approximation  determined  as 
before.  This  gave 

s0='0082, 

^  =  27,000,000, 

the  same  value  of  E  as  obtained  before.  Rougher 
observations  for  other  positions  of  the  cut  showed, 
as  one  might  expect,  that  the  s0  is  a  maximum 

when  the  cut  is  at  the  side  and  a  minimum  when  the  cut  is  at  the  top  or  bottom, 
these  positions  giving  roughly  the  same  deflection. 

*  This  small  value  as  compared  with  that  in  the  previous  experiment  is,  no  doubt,  due  to  elastic 
fatigue.  In  the  first  experiment  the  load  was  applied  by  additions  of  10  Ibs.  and  between  each  a 
considerable  time  elapsed  while  the  measurejnonts  were  being  made.  In  the  second  experiment  the 
43'7  Ibs.  load  was  put  on  and  off  at  short  intervals. 


FIG.  9. 
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Using  the  value  thus  obtained  for  E  we  get  as  our  theoretical  value  for  rl, 
when  W  =  437  Ibs.,  rZ 


Thus  the  theoretical  values  of  the  torsion  are  more  than  double  the  observed.  This 
may  perhaps  be  accounted  for  by  the  roughness  of  the  observations,  especially  of 

the  thickness  of  the  pipe,  for  it  is  evident  that,  in  the  neighbourhood  of  a0='9a, 

/Wa 
the  value  of  rl  f,T   will  change  with  extreme  rapidity.     For  example,  if  aa  =  '4:25, 

rl  --ft',.  =  214'1G  +  214'18ij  —  26770  for  77  =  5  and  the  torsions  are  then  those  in  the 

last  column  of  Table  VI  on  p.  38  agreeing  almost  exactly  with  the  observed.  Moreover 
a  tube  of  this  form  made  out  of  rolled  plate  is  probably  fairly  aeolotropic,  and  this, 
as  will  be  found,  changes  the  value  of  T  very  largely.  There  is  further  the  difficulty 
of  arranging  the  load  so  that  it  acts  through  the  centre  of  gravity,  and  above  all  the 
essential  difference  of  the  methods  of  fixing  in  practice  and  in  theory. 

§  14.     The  nature  of  the  fixing  at  the  end  of  the  beam. 

The  experimental  results  of  the  preceding  paragraph  lead  us  to  consider  more 
particularly  the  nature  of  the  end-conditions  —  what  parts  can  be  fixed  and  how,  so  as 
most  nearly  to  resemble  the  actual  conditions  of  practice.  Physically  we  see  that  any 
extended  fixing  is  impossible  with  our  solution,  for,  as  shown  in  §  9,  it  is  at  the  fixed 
end  that  the  lateral  displacements  are  greatest  and  to  them  must  be  superadded 
the  longitudinal  distortion  which  is  the  same  for  every  cross-section  of  the  beam. 
Mathematically  we  see  the  impossibility  of  extended  fixing  when  we  consider  the 
fewness  of  the  arbitrary  constants  which  arise  in  our  solution  and  which  we  alone 
have  at  our  disposal  to  describe  the  manner  of  fixing.  Bearing  our  limitations  in 
mind,  we  proceed  to  treat  some  particular  cases  of  theoretical  fixings. 

From  the  expressions  for  the  displacements  [§  4,  Equation  (20)]  we  derive  for  the 
end  cross-section,  where  z  =  0, 

W 


u=  : 


w 


.(52), 


W 


and,  arranging  the  differential  coefficients  in  a  Jacobian  scheme, 
du    du    du  W    j  W    ,  W 

dx'  dy'  dz 

v     dv     dv  W    lni  W    ,  W 

v'  di/'  dz 


dw    dw    dw 
dx'  dy'  dz 


TX  — 


W 

if 
~E~IJ  ~ 


dy 


*w 

--rrfxy  +  a.", 


.(53). 


By  means  of  the  equations  in  (52)  and  (53)  we  can  fix  certain  points  and  certain 
differential  coefficients  at  these  or  at  other  points. 
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(i)     The  fixing  at  one  point  only,  in  particular  the  centroid  of  the  end  cross-section. 

Suppose  we  aim  at  fixing  as  far  as  possible  the  neighbourhood  of  the  centroid* 
of  the  end  cross-section,  namely  (0,  y,  0).  Firstly,  to  fix  the  centroid  itself,  we 
have  from  (52), 

u=  -- 


.(54), 


, .  9w      3-y      div  .    ,  1-1-1 

and  automatically  ~— =  —=-—  =  0  tor  any  point  on  the  plane  x  —  0,   which   is  the 

(i'T  fi'JJ  fj'7 

neutral  plane. 

To  determine  a  we  may  assume  either  —  or  ^  -  to  be  zero,  since  we  know  that 

Sz         By 

a  =  a"  for  our  case  of  symmetry  about  the  axis  of  y.  This  makes  a  =  a"  =  0,  and 
means  that  there  is  to  be  no  rotational  displacement  of  the  beam  as  a  whole  about 
the  axis  of  x. 

To  determine  y  we  may  take  ;-=^-  =  0,  making  y=  —  -nTrilii.     This   really 

dy      Sx  El  '  J  3 

determines  that  there  is  to  be  no  rotational  displacement  of  the  beam  as  a  whole 
about  its  longitudinal  centroidal  line. 

From  equations  (54)  it  follows  immediately  that  a!  =  ft'  =  y  —  0. 

The  only  remaining  constants  are  ft'  and  ft",  between  which  there  is  a  relation 
given  by  either  of  equations  (37)  and  (38),  so  that  fixing  the  one  automatically  fixes' 

the  other.     We  may  either  make  —  —  0  or  7—  =  0. 

J  dz  dz 

1°.      —  =  0.     Apart  from  the  mechanical  difficulty  of  introducing  a  means  of 

fixing  which  will  fix  an  element  in  the  direction  of  z,  this  would  make  the  linear 
term  (saz)  in  the  equation  of  the  curve  of  flexure  vanish,  and  this  is  contrary  to  our 
experiments.  We  shall  therefore  make  use  of  the  second  alternative  which  was  that 
adopted  by  Saint- Venant. 

%\  ..KS  -«"_ 


2o  _      .  . 

2'     to" 

This  is  the  same  as 


since  -^  =  ^-   on   the   radius   of  symmetry,   or,   substituting    our  value   of  x   from 


equation  (30), 

r  »   /  R   \  -r 

2     Amrm'1  +  -~Am  cos  md-(&-  ft")  cos  0  +  —  -     -2  cos  26 
Lt=o  \  f     I  2  c°s  y 


*  No  property  peculiar  to  the  centroid  being  used,  the  equations  which  follow  hold  good  for  any 
point  whose  coordinates  are  (0,  y,  0)  on  the  central  axis  of  symmetry. 
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or 


giving 


Collecting  our  assumptions  and  results,  we  have  : 

If  the  centroid  of  the  end  cross-section  be  fixed  so  that 

u  =  v  =  w  =  0 


du  _  dv      div 
8x      fly      dz 


du      dv  div 

^T"  —  -T~  —  "i  ~^~  =  " 

dy      dx  dz 


•(56), 


dv  _  dw 
dz      dy 

then  the  displacement  at  any  point  (x,  y,  2)  in  the  beam  is  given  by 
W 


IF 


iv  =  ' 


...(57). 


If  we  consider  in  particular  the  curve  formed  by  the  line  of  centroids,  that  is  the 
curve  obtained  by  putting  x  =  0,  y  =  y  in  the  equations  of  (57),  we  obtain  after  a 
little  reduction 

W 


U(ll)  = 


m 


i  El  '  a  \cos  y 

assuming  uni-constant  isotropy, 
W 


8      8cosfy 


...(58). 


where  zx5  is  the  vertical  shearing  stress  at  the  centroid  of  any  cross- 
section  (from  equation  45),  and 

v(y)  =  0,         w(y)  =  0. 
(a)     Case  of  the  complete  sector. 

For  this  case  the  distribution*  of  the  vertical  shear  along  the  radius  of  symmetry 
has  been  worked  out,  and  we  can  determine  the  values  of  zx  at  the  centroid  or  at  any 

*  See  §  12,  Table  III  and  Fig.  4. 

6—2 
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other  point  with  sufficient  accuracy  by  interpolation.     We  thus  have  the  following 
table  for  the  value  of  the  coefficient  of  the  linear  term  in  u  (T/). 

TABLE  VIII. 


Angle  of  sector  : 

5° 

12° 

30° 

45° 

60° 

72° 

90° 

120° 

180° 

270° 

360° 

ZXy       Wtf 

~JT~~~EI* 

•00115 

•00674 

•04144 

•08896 

•15127 

•20726 

•30371 

•49983 

•89385 

1-20921 

oo 

It  might  be  anticipated  that  a  better  approximation  to  what  would  arise  in 
practical  fixings,  when  the  sector  has  a  small  angle,  would  be  obtained  by  taking 
the  theoretical  point  of  fixing  to  lie  between  the  centroid  and  the  circumference 


zx 


end  of  the  radius  of  symmetry.     The  nature  of  the  corresponding  values  of    -  can  be 

immediately  seen  from  Table  III.     Such  a  point  of  fixing  would  enable  us  to  avoid 
the  infinite  linear  term  in  the  deflection  for  the  split  cylinder. 

(I)}     Case  of  the  sector,  complete  or  curtate,  when  y  =  Zir. 

In  this  case  y  is  at  the  origin  and  when  y=0,  zXy  becomes  infinite  for  the  com- 
plete sector.  Moreover  we  must  always  have  the  relation  of  (38)  satisfied  whatever 
be  our  theoretical  method  of  fixing  and  it  is  clear  that  for  our  case  of  y  —  2vr  the  value 
of  /3  —  ft"  will  have  to  be  infinite,  since  r  is  not  zero. 

These  mathematical  infinities  need  not  deter  us  from  applying  the  theory  to 
physical  problems.  In  the  actual  experiments  we  cannot  get  angles  of  exactly  2;r  and 
we  have  no  means  of  practical  fixing  in  accordance  with  the  theoretical  methods.  In 
fact  the  true  interpretation  of  these  infinities  is  that  we  need  not  be  surprised  at 
relatively  large  numbers  arising  at  the  corresponding  points  in  physical  problems  and 
it  was  the  large  value  found  in  the  experiments  with  the  cut  pipe  which  led  us 
to  the  investigations  in  this  section. 

(ii)     fixing  at  two  points — case  of  cylindrical  split  pipe. 

Instead  of  fixing  one  point  and  certain  directions  through  that  point  we  may  use 
the  arbitrary  constants  in  the  solution  to  fix  two  points  of  the  terminal  face.  The 
choice  of  the  two  points  would  naturally  vary  with  different  kinds  of  sector  and  we 
shall  devote  our  attention  immediately  to  the  cylindrical  split  pipe. 

We  shall  choose  the  two  points  ( +  a,  0,  0),  i.e.  the  two  points  where  the  plane  of 
flexure  crosses  the  outer  circle  of  the  end-section.  The  equations  for  u  =  v  =  w  =  0 
will  then  be 


W 


iv  =  • 


r  =  a 


X(0  = 


+  a'=0 
±ya   +  /3'  =  0 

+  (3"a  +  y  =  0 


.(59). 
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From  these  we  deduce  immediately 

W  , 


and 


To  determine  a  we  may  put  ^-  =  0  at  (±a.  0,  0) — making  a  =  0  and  virtually 
preventing  a  rotation  about  the  diameter  through  (±ci,  0,  0). 
The  equation  for  ft"  is  from  equation  (30) 

IT7 

r[i  +  (i-1 


a  = 


+         sin         -  (ft  -  ft     a  - 


or, 


whence  ft  =  1 

or,  substituting  for  Am  and  Bm  from  (34), 


.(60), 


P  ~ 


2ir 


8x520-6(1+1?) 
m  (m2  — 4) 


m(m2-9) 


7T 


a(n(2i+l)- 


With  these  determinations  of  the  arbitrary  constants  the  displacement  equations 

become  ^ 

u=  - ryz  +        [l 77  (a;2 - 1/2)  (l 


W 


...(62), 


iv  = 


- 


and  in  particular  the  curve  of  centroids,  which  for  the  sector  in  which  y  —  2ir  is  the 
same  as  the  curve  of  centres,  is 


W  W 

—      -  fi/z2  —  x  Z3~l  -4-  Rz  —  — 
"~  EI  6     -^ 


.(63). 


The  last  constant  term  is  of  course  the  result  of  the  anticlastic  curvature  and  as  we 
shall  show  is  of  no  importance  compared  with  the  value  of  ftl.  The  value  of  ft  is 
identical  with  s0  as  in  all  cases  where  the  centroid  coincides  with  the  centre. 

Returning,  therefore,  to  our  value  of  ft  in  equation  (61)  we  find  that  the  first 
ten  terms  of  the  series  are  needful  to  give  adequate  accuracy  when  a0  =  ^a,  but  the 
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only  terms  that  are  of  real  importance  are  those  containing  the  large  constant 
520'6  (1  +  17).  The  result  of  the  calculation  is  that 

927-931  ...........................  (64), 

an  extraordinary  value  when  compared  with  the  values  (for  angles  other  than  360°) 
in  Table  VIII,  and  with  the  comparable  values  for  the  case  of  the  beam  of  rectangular 
cross-section  calculated  by  Saint-  Venant*. 

This  large  value  is,  however,  confirmed  by  our  experiment  described  in  §  13.  For 
when  we  calculate  ft  for  a  pipe  of  the  dimensions  used  in  the  experiment,  we  find 
that  /3  =  '0240.  This  is  somewhat  larger  than  the  observed  value  which  was  '014. 
However  it  is  of  the  same  order  when  we  compare  it  with  other  cases. 

It  is  interesting,  for  instance,  to  compare  it  with  the  value  of  ft  deduced  by  the 
same  method  of  fixing  if  the  pipe  be  closed.  It  is  easy  to  show  for  this  case  of  the 
complete  hollow  circular  cylindrical  pipe  that 

W  I  n^n  2\  W 

6-(ft-ft")rsme  ...(65), 


the  notation  being  the  same  as  used  for  the  curtate  cylindrical  sector,  and  that  the 
value  of  ft,  when  the  points  r  =  a,  6=  ±H>  2  =  0  are  fixed,  is  given  by 

W 

K)  ..............................  (66), 


or,  when  a  =  0'50  inches   and    c»0  =  0'45  inches, 

£=•000053, 
—  about  g-^-yth  of  the  corresponding  ft  when  there  is  a  cut  in  the  pipe. 

Now  the  plugging  of  a  split  pipe  and  the  clamping  of  the  plugged  end  in  the 
chuck  of  a  lathe  may  not  very  closely  correspond  to  fixing  indefinitely  small  elements  at 
the  top  and  bottom  of  the  external  circumference  at  the  fixed  end  of  the  pipe.  It 
is  worth  while  inquiring  what  value  ft  would  take,  if  we  fixed  external  elements  on 
the  horizontal  diameter  at  the  clamped  end.  We  may  take  as  our  conditions  : 

u,  v,  w  =  0  for  z  =  x  =  0,  r  =  a,  6  =  0, 
and  again  for  r  =  a,  6  =  IT. 

W 

From  (35)  we  have  u  —  0  =  ^  (  —  ^r?a2)  I  ±  ya  +  a', 


whence  flow  y  =  0,  a'  =  ^    „     r)l,  v  =  0,  ft'  =  0  ; 


—  0  =  x  (a,  IT)  —  a"a  +  y. 


*  "Memoir  on  Flexure,"  §  29  (Journal  de  Lionmlle,,  Series  II,  Tome  i  (1856)).    Quoted  in  Todhunter 
and  Pearson's  History,  Vol.  n,  Section  [96]. 
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Now  x(a>  0)  'ls  zero. 


»=o 

does  not  involve  any  additional  unknown  constant. 
Hence  these  equations  determine 

a"  =  2a  X  (a>  *}  and  7'  =  ~  ^X  (°"  ^' 

Tt  still  remains  to  select  a  condition  which  provides  /3.     A  suitable  limitation  seems 
to  be  that  dw/dx  should  be  zero  at  non-split  edge,  or 

_dw_(dx\  W 

--=a~~ 


00  /  B   \  Was 

which  leads  to  ft  =  2    Amam~l  +    „  ^  )  m  +  TO,  -  1  17    ,rr  . 

i=o\  a      I  rLil 

Wa 

Substituting  for  Am  and  Bm  from  (34)  and  the  value   ^  x  520'6  (1  +  rj)  for  T,  we  have 


I, 


(-I)'    -8x520-6(1 

Hi3 -4 


3- 


Wa? 

or,  if  T,  =  I,  /8  =          x  959-675, 


a  result  comparing  fairly  with  the  previous  high  value  927'931  obtained  by  fixing 
elements  on  the  vertical  diameter.  It  seems  necessary  to  infer  that  with  perfect 
clamping  we  should  reach  a  value  of  /8  comparable  with  that  given  in  (64). 

Experimental  verifications  of  the  foregoing  theory. 

It  will  be  of  advantage  at  this  stage  to  summarise  some  of  the  results  of  experi- 
ments made  with  a  view  to  increasing  the  faith  in  the  theory.  In  all  the  cases 
considered  here,  the  deflection  (Sz)  was  observed  at  various  points  (z)  along  the  lowest 
generator*  of  the  cylinder,  and  the  constants  s0  and  s1  in  the  formula 

8,  =  stz  +  Sl(ite-W  ..............................  (67) 

were  deduced  from  the  measurements  by  the  method  of  least  squares.  The  results 
are  given  in  Table  IX.  Remembering  that  s0  seems  naturally  to  arise  in  the  form 

s0  =  (  Wd'/EI)  x  a  quantity  depending  on  the  shape  of  the  cross-section, 

Wa2 
*  sa  will  therefore  differ  from  /J  by  a  term  of  the  form  -=j  £77.     This  term  would  have  no  sensible 

effect  in  the  case  of  the  split  pipes,  but  would  mean  that  s0  is  less  than  /?  by  about  C  °/o  in  the  case  of  the 
non-split  tube. 
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where  a  is  half  the  typical  vertical  dimension  of  the  cross-section  of  the  beam*,  we 
shall  write,  after  Saint- Tenant,  the  total  deflection  as 

(68). 


The  quantity  cr0  is  tabulated  as  well  as  the  ratio  sa/sl.     For  purposes  of  com- 
parison, of  course,  <TO  is  superior  to  s0/s1. 

TABLE  IX.    Experimental  results  on  the  nature  of  the  deflection 
with  various  beams :  8Z  —  saz  +  s1  (^Iz*  —  ^'/). 


Description  and  dimensions  of  beam 

»u 
(W 

*i 

£0 

Si 

«0/*1 

°~  a" 

Stretch-modulus 
E  Ibs.  per  sq.  in. 

Split  pipe,  (soft  steel)  : 

a  =0-50  in. 

cut  at  side 

»„  =  0'45  in. 

•0136 

•000096 

140 

560 

27,000,000 

I  =  15  in. 

(-0240) 

Same  pipe  :  cut  on  top 

a  =  0-50  in. 

a0  =  0'45in.     i  -0082 

•000096 

85 

340 

27,000,000 

I  =15  in.         (-0068) 

Soft  steel  pipe  without 

a  =|  in. 

cut 

«o  =  TF  in- 

•00181 

•000236           7-7 

54 

23,000,000 

I  =9-5  in.        (-000061) 

Bass-wood  cylinder 
with  saw-cut  to 
central  axis 

© 

a  =  |  in. 
I  =32  in. 

•0151 
(-00046) 

•0000681 

220 

570 

1,470,000 

The  experiments  were  only  rough  and  little  was  expected  of  them  beyond  some 
slight  verification  of  the  theory.  The  values,  however,  of  the  Young's  Modulus,  E, 
are  in  accordance  with  what  we  should  anticipate  for  the  materials  used.  The  chief 
divergence  lies  in  the  values  of  s0  when  compared  with  the  theoretical  values  of  /8  in 
the  case  of  the  solid  pipe  and  the  wooden  cylinder  with  the  saw-cut.  But  the  nature 
of  clamping,  by  a  chuck  for  the  uncut  pipe  and  a  simple  clamp  for  the  bass-wood 
cylinder,  is  really  not  that  of  two  fixed  points  and  may  well  provide  somewhat  more 
closely  the  high  values  of  a  central  fixing :  see  (b),  p.  44. 

We  may  note  :  1°.  the  difference  between  the  tr0's  for  the  split  pipe  when  the 
cut  is  at  the  side  and  when  it  is  on  the  top  (the  arrows  in  the  drawings  of  the  cross- 
sections  being  in  the  direction  of  the  flexure) ; 

2°.  the  small  value  of  <TO  for  the  complete  tube  although  it  is  not  so  small  as  the 
method  of  fixing  by  the  two  points  ( +  a,  0,  0)  described  above  would  provide  ; 

3°.  the  cr0  for  the  wooden  beam  with  the  saw-cut  up  to  the  central  axis.  This 
is  the  physical  approximation  to  the  non-curtate  sector  of  angle  2-rr.  The  o-0  is  large 
as  we  should  expect  from  a  consideration  of  the  order  of  the  ratio  of  stretch-modulus 
(E}  to  slide-modulus  (ju.)  for  wood.  This  is  considered  more  fully  in  the  next  section  of 
the  paper. 

*  In  our  case  the  outer  radius  of  the  sector,  in  the  case  of  a  beam  of  rectangular  section  the  semi- 
thickness,  and  for  an  elliptic  section  the  semi-axis  parallel  to  the  direction  of  flexure. 
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IV.     SOLUTION  FOR  CERTAIN  AEOLOTROPIC  CASES. 

§15.     A  solution  for  aeolotropic  material. 

Reverting  to  our  original  equations  in  §  3  for  a  uniformly  aeolotropic  substance, 
we  see  that  the  main  change  in  the  analysis  that  a  consideration  of  aeolotropy 
requires  will  be  introduced  by  the  equation  (15)  of  p.  10 


instead  of  the  usual  Laplace  equation.  It  is  natural  to  try  a  projection,  which  will 
transform  this  equation  into  Laplace's  equation.  Let  us  project,  then,  with  the 
substitution 


This  projection  depends  on  the  nature  of  the  aeolotropy,  circles  being  changed 
into  ellipses  of  semi-axes  proportional  to  J^  and  ,//*„,  but  in  such  a  way  that  there 
is  no  change  in  the  lateral  dimension  of  the  beam. 

If  the  boundary  be  <j>  (x,  y)  =  0, 

the  transformed  boundary  will  be 


or  p  (x',y')  =  (j, 

if  we  agree  to  denote  everything  in  connection  with  the  transformed  problem  with 
dashed  letters.     Also 

dx' 
cos  (x1 i/)  = 


+ 
-^/=«.(»)x    /.^.    ™v (70,, 

and,  similarly, 

COS(,V')=OOS(,,,)X 

I 


Thus  the  problem  of  the  aeolotropic  case,  namely : 
To  find  x  (x,  y]  and  T  to  satisfy 
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with  the  boundary  condition 

u, .--  cos  (xv)  +  /A2  ~  cos  (i/v)  =  r(fry  cos  (xv)  —  frx  cos  (yv)) 

ox  oy 

W 

„         ,„          .„_,         _.,,,,„,,  ^,      -<  CQ8lxv\  +  IJi;_7)    „    \a;y  COS 


-/A!  (/3-/3")  COS  (ajj^+^^-a")  cos  (yv)     (72), 

and  the  Conditions  I  and  II, 

0=  \\oo8(yv}d* (73), 


is  projected  into  the  problem  : 

To  find  x'  (d,  y')  and  T'  to  satisfy 


>j 

with  the  boundary  condition 

r^  cos  (xV)  +  ~j  cos  (yV)  =  /  /—  (y'  cos  (a;'^)  —  x'  cos  (y V)) 
oic  oy  \j  p.., 

W      ILL,  r/i?iMi     i«      £* — T?!/*!  —  2rj.,/A.,    ,A          /    ,    ,,.       -fc  —  ^lA1]     /    /         /   /  /\  I 

+  TTT    /—         -«2  +  -  — — "  w     cos  (*V )  +  —         -aryooBiyy) 

jfcj  \j  fr  |_\  2/t2  2tl2  /  LL, 

-((3-(B")f     ^cos(x?v')  +  (a.-a.")'cos(y'v')   (75), 

\/  fr 

and  the  Conditions  I  and  II, 

n         I      /  t    I    f\  J  I  /'7C\ 

0  —   I  V    COS  I  II  V  I  US ('O)> 

i  A  V^/        /  \        / ' 


+    /A 

-         (E  -  ;!  n  a  ^-'  x'27/  +  -          -..    /3 

t        tW«   ^+'  2/A2 

^dx'dy*     ...............  (77). 


This  projected  problem  is  analytically  exactly  the  same  as  for  the  isotropic  case*. 
The  only  difference  from  the  latter  lies  in  the  different  constants  in  the  three  condition- 
equations.  Hence  any  solution  that  has  been  obtained  for  an  isotropic  case  can  be 
used  to  obtain  a  solution  for  a  connected  problem  in  an  aeolotropic  material.  There 
is  however  the  limitation  which  must  be  emphasised1  (to  the  depreciation  of  the 
method),  that  the  form  of  the  soluble  cases  depends  essentially  on  the  solubility 


*  Note  that  /  =  (i^lit^ft  /'  has  been  retained  and  is  therefore  =  i  (f^//^)*  (a4  -  a04)  (y'  -  sin  y')  for  the 
'  curtate  elliptical  sector.' 
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of  a  corresponding  isotropic  case  and  also  on  the  nature  of  the  aeolotropy  —  to  be 
precise,  on  the  ratio  of  the  two  slide-moduli  ^,,  /z2.  The  result  is  that  the  soluble 
cases  cannot  always  be  for  cross-sections  of  common  occurrence,  but  must  be  for 
distortions  of  such.  In  the  uncertainty  of  the  elastic  constants  of  aeolotropy,  how- 
ever, this  disadvantage  may  in  many  cases  be  of  little  importance. 

§  16.      Case  of  the  curtate  sector  cross-section. 

Let   us  take  the  projected  figure  to  be   the  curtate  circular  sector  which    is 
enclosed  between  the  circles 


(78). 


and  the  lines  yf=  +  «'tan  \y 

Then  on  projecting  back  by  the  transformation 


the  cross-section  given  will  be  a  '  curtate  elliptical  sector  '  enclosed  between  the  ellipses 


af +«•««' 

Mi 

— 


.(79). 


and  the  lines  y  =  ±      —x  tan  ly'  =  ± x  tan  \y,  say 

V  Mi 

Following  the  method  of  treatment  used  for  the  isotropic  case  in  §  6,  we  get  the 
following  solution  for  ^'  in  the  projected  problem  : 


x'  =  i  I A'mr'm  +  5^1  sin  mff -(ft-  ft"}'    fa  r' sin  6'  +  r'    fa  ^-  -.  sin  20' 

»=o  \  r  /  v  M"  v  M"  2  cos  y 

[ 


where 


1  W_     Mrt^-i'h/*.-: 

3EI*J  p.. 


COS  I/ 


/i.  / 

4'     =-«T'       /&      —  lZ 

V  /*,  (m«  - 


™m+2 

~ 


m+2 


4)  my'    a2™  - 


1008|y//f^-»1u,       1         .4  J?- 971,0,      1     \am+3-a( 


/„  /  _  1  V         am+'*ri  m+2  (nm~* 

z>/     _  o   /      /Mi         V       J  )  ao         l.a 

\/  fL,  (ma  -  4)  my'  a2™  -  a,,2™ 


m2—  1 

,  m-2\ 


my 


m2-! 


m2-9 


„«!+:)„  m+3  /.-,Mj-3         „  m-3\ 
01 ao          (<l          ~Oo / 

„•««        ,v  'J;/i 
t*       —  (*„ 


and 


, 


.(81). 
7—2 
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Condition  I  will  as  before  (p.  17)  be  identically  satisfied  for  symmetry  perpen- 
dicular to  the  plane  of  loading  applying  Condition  II,  and  remembering  the  value 
of  /  (see  ftn.  p.  50),  we  get  as  equation  for  the  torsion  T', 

r'  [(a4  -  a/)  (yf  -  tan  y'}  +  I   ,—  —, ( —.    -^^ 

,-=0  (m~  —  4)  (m  +  2)  my      eP*  —  «„ '" 

oo  CA  n*n*lnm~- n  m~- 

„  o4  a  (MI*  \a       — ua 


,Co  (m*  -  4 )  (m  -  2)  my' 


W  T  5      2cos-iy'    /QE— An.it.       1  (i^-SBn.u,      1 

V  ?Lt  /-r~-  /-r~- 

—  TTf       * 


i =o  ( m  +  2)  my'  V        /LI,         m-  -  I  /x,          m2  - 


X 


•      2  cos 


X 

2   E  —  6r7.,u., 

'-^2  cos  iv'  tan  3  y  (a5  -  ans) 


v/i-'  -  1               /*„           m2  -  9 
r/5r/  5  ( sim~'A n  m~ 

\,V    Lt-Q      I  \-v  L*Q 


»'-<) 

...........................  (82). 


Now  projecting  back,  so  far  as  this  equation  is  concerned,  consists  merely  in 
replacing*  r'  by  T,  and  y'  by  y,  where  y'  =  2  tan"1  (    /—  tan  |y)  and  retaining  /,  so 

\  \.  '    LL  ,  I 


vv 

that  with  these  changes  the  equation  (82)  gives  the  value  of  the  torsion  accom- 
panying the  flexure  of  a  beam  having  the  curtate  elliptical  sector  cross-section  given 
by  the  equations  (79). 

§  17.     Solution  in  the  case  of  transverse  isotropy. 

Transverse  isotropy  consists  in  the  slide-modulus  and  the  Poisson  ratio  being 
independent  of  the  direction  perpendicular  to  the  longitudinal  axis  of  the  beam.  We 
may  write  /HJ  =  JU,,  =  /A  and  Yi1  =  7).,  =  r),  but  we  do  not  necessarily  have  .Z?  =  2/n  ( 1  +  77) 
which  would  give  the  case  of  complete  isotropy. 

Applying  the  method,  of  the  previous  section  we  see  that  the  projected  figure  is 
the  same  as  the  original  figure  in  this  case,  but  we  can  deduce  the  result  from  the 
aeolotropic  case  and  have  immediately  from  equation  (82), 


and   I. 


*  Of  course  'replacing'  T'  by  r,  and   retaining  7  does   not   involve   equality  of   T'  and  T  or  of   /' 

r 
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/  4  64  (am+2-a0m+2)2 

r    (a  -  <)  (y  -  tan  y)  +  -  -  - 


»  64 


."" 


i=o  (w1  -  4)  (wi-  2)  my      ~  a""  -a 

JF  T  •      2  cos  //    £          \       1  IE        ,  \       1 


A'/          '    1 1     _i_  >,>\  \\  *  I       -       1         \  II       a       Q 

/,.m  +  2        „  »»  +  a\  I r.m+3        „  m  +  3\ 

fa      -  «„      j  (a      -  a0      ) 


6  --  47i  —    6  -- 


, 
m-—l      \    p.  '/nr  —  9 

a5a0°  (a"1"''  -  a,,"1-3)  (am~2  -  <"~2) 


(83). 


y 

A  comparison  with  (82)  shows  us  that  this  result  is  absolutely  the  same  if  the 
appropriate  JS/fJ^  and  rj.2  be  used.  In  other  words  the  difference  of  the  two  cases  lies 
solely  in  the  distortion  of  the  auxiliary  cross-  section. 

§  18.     Note  on  the  torsion  of  an  aeolotropic  beam  ivith  flexure.     Of.  §  9. 
The  equation  for  the  vertical  displacement  at  the  point  (x,  y,  z)  in  the  aeolotropic 
case  is  (equation  14) 


the  only  change  from  the  isotropic  case  being  the  introduction  of  the  two  Poisson's 
ratios  T/,  and  17.,.  When  we  consider  the  displacements  on  the  horizontal  axis  of 
symmetry  the  only  Poisson's  ratio  occurring  is  the  77,,,  and  we  have  as  the  torsion  in 
the  case  of  the  beam  discussed  in  the  previous  section 


a}  ........................  (84). 

Compare  equation  (42). 

Thus  the  torsion  due  to  the  anticlastic  curvature  depends  only  on  that  Poisson's 
ratio  which  is  the  contraction  in  the  direction  of  y  (i.e.  perpendicular  to  the  plane  of 
flexure)  due  to  an  extension  in  the  direction  of  z,  the  longitudinal  axis. 

§  19.     Applications  to  certain  aeolotropic  materials  and  experimental  results. 

The  interest  in  the  aeolotropic  case  does  not  lie  so  much  in  the  immediate 
usefulness  of  the  formulae  as  in  the  estimation  of  the  order  of  magnitude  of  the 
changes  that  the  aeolotropy  of  a  commonly  occurring  kind  is  likely  to  produce. 
Aeolotropic  materials  other  than  those  of  regular  crystalline  structure  are  by  nature 
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so  variable  in  their  aeolotropy  that  it  is  impossible  to  give  dependable  values  to  their 
elastic  constants.  In  the  case  of  the  common  woods  for  instance  all  depends  on  the 
part  of  the  tree  trunk  that  is  selected,  and  in  any  one  plank  the  aeolotropy  may  be  so 
variable  as  to  upset  even  the  nature  of  the  mathematical  result  proved  for  uniform 
aeolotropy.  Thus  in  an  experiment  that  was  made  with  a  plank  about  four  inches 
broad  and  bevelled  so  as  to  have  for  cross-section  approximately  a  circular  sector  of 
small  angle,  it  was  found  that  the  torsion  caused  by  flexure  was  most  pronounced  but 
in  the  direction  contrary  to  that  given  by  theory.  The  annual  rings  in  the  wood  were 
noticed  to  be  much  wider  apart  towards  one  edge  of  the  plank  and,  this  explanation 
having  suggested  itself,  the  experiment  was  repeated  with  the  'complementary'  piece 
of  wood,  i.e.  with  a  plank  in  which  the  rings  broadened  towards  the  other  edge  of  the 
bevelling.  This  time  the  torsion  was  of  the  same  sign  as  in  the  mathematical  theory 
but  the  •  magnitude  was  quite  out  of  proportion  to  what  might  be  expected  with 
uniform  aeolotropy. 

It  is  conceivable  that  this  property  of  the  variability  of  the  -elastic  constants  in 
wood  might  be  made  use  of  in  those  cases  of  practice  where  it  is  important  that  the 
flexure  of  a  beam  of  asymmetrical  cross-section  should  cause  little  or  no  torsion  or  a 
torsion  in  a  certain  direction.  Although  it  might  be  difficult  to  select  single  planks 
that  would  give  the  suitable  torsion,  it  might  be  easier  to  make  compound  torsion- 
free  beams  of  several  layers  of  different  pieces  of  wood: 

As  a  rough  guide  to  the  magnitude  of  the  change  due  to  aeolotropy,  it  is  easy 

,  E  . 
to  consider  the  case  of  transverse  isotropy  (§  17).     If  --  is  fairly  large  so  that  we  can 

P- 

Tjl 

neglect  multiples  of  17  in  comparison  with  — ,  we  see  from  equation  (83)  that  the 

J?  T?  f* 

value  of  T  is  roughly  proportional  to  —  .    The  value  of  --  for  an  isotropic  material  is  -  , 

P  /* 

but  for  wood  it  is  much  higher,  if  we  take  E  in  the  direction  of  the  wood  grain.     In 

E 
the  case  of  beechwood,  for  example,  we  can  estimate*  that  -  -  is  about  6,  and  for 

j* 

pinewood  about  30. 

Some  measurements  were  made  of  the  torsion  produced  in  a  round  beam  of  bass- 
wood  having  a  saw-cut  through  to  the  central  axis.  The  beam  had  then  approxi- 
mately the  cross-section  of  a  complete  circular  sector  of  angle  360°.  The  radius  of 
the  beam  was  |-  of  an  inch  and  the  length  of  the  flexed  part  about  36  inches.  The 
flexing  weight  which  was  used  varied  up  to  12  Ibs. 

To  measure  the  amount  of  torsion  one  end  of  a  stiff  wire  was  inserted  in  the  saw- 
cut  at  the  free  end  of  the  beam  so  that  the  wire  stuck  out  horizontally  and  at  right 

*  The  estimation  is  based  on  the  values  of  the  stretch-modulus  given  by  Wertheim  and  Chevandier 
(Comptes  Rendus,  T.  xxm,  p.  663,  1846)  and  of  the  Poisson's  ratios  found  by  Mallock  (Proc.  Roy.  Soc., 
Vol.  xxix,  p.  157,  1879),  and  is  made  on  the  assumption  of  the  'partial  ellipsoidal'  elasticity  of 
Saint- Venant  (cf.  Todhunter  and  Pearson's  History,  Vol.  II,  Part  I,  §  [314]). 
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angles  to  the  beam's  axis  when  the  beam  was  in  its  normal  position.  The  rotation  of 
the  wire  pointer  was  measured  by  noting  the  change  in  the  position  of  its  shadow  cast 
on  a  piece  of  sectional  paper  by  a  light  at  some  distance  behind  the  fixed  end  of  the 
beam. 

Without  going  into  the  details  of  the  measurements,  it  was  found  that  the  total 
torsion  amounted  to  '025  radians  (or  1'5  degrees)  with  a  flexing  weight  of  12  Ibs.  and 
a  deflection  of  about  1^  inches. 

For  an  aeolotropic  split  cylinder  the  accurate  torsional  relation  is 


,-S(4*+%)x -10811, 


,58733. 


In  the  present  case  T=  '025/36,  E=  1  '470,000,  a  =  |  and  /=  '11984,22489.  Thus 
we  have  158 '189(5  =  ^E/^  +  r).,.  The  annular  rings  were  approximately  parallel  to  the 
plane  of  flexure,  and  if  we  assume  77,  =  '45,  we  have  E/p.^  =  39 '44,  a  high  value,  but  not 
without  parallel  in  the  case  of  other  woods. 

For  most  practical  purposes  this  method  of  transverse  isotropy  is  as  exact  as  our 
knowledge  and  the  variability  of  the  elastic  constants  of  wood  entitle  us  to  use,  but 
it  is  worth  while  to  investigate  one  instance  more  closely  with  a  more  complete  system 
of  elastic  constants,  and  to  see  the  differences  which  arise  according  to  the  plane  of 
aeolotropic  symmetry  which  is  taken  as  the  plane  of  flexure. 

Using  the  figures  of  Wertheim  and  Chevandier  and  of  Mallock  the  following 
systems  for  beech  wood  and  pine  wood  were  arrived  at  on  the  assumption  of  the  'partial 
ellipsoidal'  elasticity  of  Saint- Venant. 

TABLE  X.     Elastic  Constants  of  Beechwood  and  Pinewood  on 
the  basis  of  '  partial  ellipsoidal '  elasticity  *. 


Elastic  Constants 

Beechwood 

Pinewood 

E 

980 

1113 

•t]  (perpendicular  to  annual  rings) 
rj'  (parallel  to  annual  rings) 

•408 
•530 
172-9 

•372 
•486 
45-3 

P-" 

132-7 

27-1 

n  :  the  coefficient  of  '  partial  \ 
ellipsoidal  '  elasticity        / 

•754 

•150 

Of  course  no  great  reliance  can  be  laid  on  the  accuracy  of  these  figures,  but  they 
provide  systems  such  as  are  likely  to  occur  in  wood  and  give  some  idea  of  how  in 
practice  the  aeolotropic  solution  given  above  in  §  16  would  work  out. 

The  case  dealt  with  for  the  isotropic  solution  at  y  =  12°  was  selected  for 
investigation  aeolotropically — that  is,  the  'projected'  problem  was  taken  to  be  that 

*  Moduli  measured  in  kilogrammes  per  sq.  millimetre. 
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concerning  a  non-curtate  circular  sector  of  12°  angle.     The  results  of  the  calculations 
are  given  in  the  following  summary  : 

TABLE  XI.     Aeolotropie  case :  Numerical  results. 


Material 

Angle 

IWa 

•fa 

IWa  . 
^  1  —  -  for  isotropic 

/  at 

case  with  the  same 

1°     Isotropic  with  t\  —  \        ...          ...          ...          

12°0' 

•2846 

angle  as  for  the 
aeolotropic  case 

and  ?j  =  J 

2°.    Beechwood  : 
(a)    with  annual  rings  iu  plane  of  flexure 
(6)    with  annual  rings  perpendicular  to  plane  of  flexure 

13°  40'-9 
10°31'-3 

•6828 
•5253 

•2837 
•2852 

3°.    Pinewood  : 

(a)    with  annual  rings  in  plane  of  flexure 
(b)    with  annual  rings  perpendicular  to  plane  of  flexure 

15°  28'-6 
9°  17'-7 

1-9410 
1-2226 

•2825 

•2855 

For  the  purposes  of  comparison  of  these  aeolotropic  cases  with  the  isotropic  cases  of 
a  circular  sector  of  equal  angle,  it  is  clear  that  the  elliptical  shape  of  the  aeolotropic 
sectors  is  not  of  great  importance.  This  is  illustrated  in  Fig.  1 0. 


y'=  12°    0'.    Projective  solution,  true  circular  sector. 

y  =  15°  29'.     Pinewood,  annual  rings  in  plane  of  flexure,  sector  with  elliptic  arc. 

y  =    9°  18'.     Pinewood,  annual  rings  perpendicular  to  plane  of  flexure,  sector  with  elliptic  arc. 

FIG.  10.    Diagram  to  illustrate  closeness  of  elliptic  sectors  to  true  circular  sectors. 

The  most  noteworthy  results  of  these  calculations  are  to  show  (i)  the  etf'ect  of 
the  direction  of  the  annual  rings  in  the  wood  on  the  torsion  due  to  unbalanced  shear 
couple,  and  (ii)  the  great  increase  in  the  flexural  torsion  when  we  pass  from  a  material 
such  as  steel  to  wood.  They  serve  to  emphasise  the  difficulty  of  making  applications 
of  the  mathematical  theory  to  give  numerical  predictions  regarding  the  behaviour  of 
such  a  variable  substance  as  wood.  Without  more  precise  determinations  of  the  elastic 
constants  of  any  individual  plank  we  can  at  best  give  descriptive  properties  only  and 
even  then  may  be  disappointed  not  to  find  them  fully  confirmed  by  experiment ! 
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V.     SUMMARY  OF  ANALYTICAL  RESULTS  OBTAINED  FOR  BEAMS  OF  CROSS-SECTIONS 
OF  ASYMMETRICAL  FORMS  OTHER  THAN  THAT  OF  THE  CIRCULAR  SECTOR. 

In  tins  section  are  collected  the  analytical  discussions  of  the  problem  of  flexure 
for  various  other  asymmetrical  cross-sections,  which  would  appear  to  be  amenable  to 
arithmetical,  if  not  to  complete  analytical  solution.  The  discussion  is  confessedly 
incomplete  in  that  the  results  obtained  are  chiefly  in  analytical  forms  which  can  convey 
no  immediately  apparent  physical  meaning.  To  reach  the  stage  of  physical  applica- 
tion would  require  months  of  calculation  arid,  as  long  as  the  exigencies  of  more 
pressing  war-work  do  not  permit  of  the  necessary  leisure,  we  can  but  look  forward  to 
the  arithmetical  completion  of  the  work  in  a  continuation  of  the  present  memoir.  It 
seems  however  desirable  to  describe  here  what  has  been  done. 

§  20.  Case  of  a  cross-section  of  algebraic  equation  derived  from  assuming  an 
integral  algebraic  form  for  the  function  x  (%,  y}-  An  attempt  to  extend  a  method 
used  by  Saint-  Venant. 

In  looking  for  guidance  in  the  works  of  Saint-Venant  it  is  natural  to  expect 
that  soluble  cases  might  be  found  in  developments  of  his  algebraic  method  of 
solution*.  The  basis  of  the  method  is  to  start  with  the  assumption  that  the  function 
X  (x,  y),  which  arises  in  the  expression  for  the  longitudinal  displacement  w\,  is  an 
integral  algebraic  function  of  x  and  y,  and  to  find  therefrom  the  boundary  of  the  cross- 
section  for  which  this  function  gives  a  consistent  solution. 

Thus,  if  we  assume  that 


.....................  -(85), 

we  find  well-known  relations  between  the  coefficients  which  must  hold  if  the  equation 

^    1v      '       3X^=o  •' 

dx2     dy* 

is  to  be  satisfied.      When  we  substitute  the  resulting  harmonic  expression  for  ^  in 
the  boundary  equation  (22) 

xz  cos  (xv]  +  yz  cos  (yv)  =  0, 
that  is,  in  terms  of  the  differentials  dx,  dy  along  the  boundary  of  the  cross-section 

xzdy  —  yzdx  =  0     ..............................  (86), 

we  have  the  differential  equation,  of  the  boundary  of  the  cross-section. 

Saint-Venant  found  that  a  large  number  of  interesting  cases,  including  that  of 

*  See  the  Memoire  sur  la  flexion  (cited  on  p.  4),  §§  18  —  23. 
t  Cf.  equations  (14)  or  (20)  above. 
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the  elliptic  cross-section,  could  be  discussed  by  retaining  only  the  harmonic  x3  —  3xif. 
This  leads  to  boundaries  contained  in  the  general  equation* 


.....(87), 

where  m,  c,  c0  c2  are,  with  certain  restrictions,  constants  at  our  disposal.  We  cannot, 
however,  make  use  of  this  equation  for  any  but  symmetrical  cross-sections,  because 
in  Saint-  Venant's  form  of  the  shears  yz  and  xz  the  torsion  constant  r  was  taken  to  be 
zero.  Thereby  Saint-  Venant  was  virtually  assuming  symmetry  of  cross-section,  and 
he  seems  to  have  overlooked  this  consideration  in  discussing  f  the  geometrical  forms 

fyY) 

of  contours  which  would  be  given  by  odd  values  of  —    -  in  equation  (87).     He  does 

not  however  proceed  to  discuss  in  detail  the  flexure  of  a  beam  having  any  such 
asymmetrical  cross-section. 

We  have  therefore  to  follow  up  Saint-Venant's  treatment  with  the  more  general 
form  of  equation  (86)  that  is  given  in  §  4  above.  Omitting  the  details  of  the  work, 
we  have  the  result  that  if  we  take 

W 

......  (88), 


the  boundary  of  the  cross-section  which  corresponds  is  given  by 
j»_       W.  ,  1— m       W(l—m(l- 


: (89). 

The  latter  equation  may  be  written  as 


m 

(,,\l^m        ™2  ..        .f. 

f)  +l-ff4=°  ...........................  <90> 


where  m,  £,  £J,  a,  b  are  constants  which  are  not  at  our  disposal  until  Conditions  I 
and  II  of  §  4  have  been  satisfied.  Condition  I  is  easily  shown  to  be  immediately 
satisfied,  but  Condition  II  presents  greater  difficulties,  for  the  boundary  as  well  as 
the  integrand  of  the  surface  integrals  involved  depend  on  the  constants  of  equation 
(90).  As  that  equation  stands,  x  is  expressible  in  terms  of  a  trinomial  in  y  and  it 
appears  that  in  order  to  render  the  equation  workable  we  must  assume  one  of  the 
coefficients  to  be  small  and  so  be  able  to  expand  in  powers  of  this  coefficient,  the 
integrals  of  Condition  II  then  being  reducible  to  B-functions.  This  was  attempted 
in  two  cases,  (i)  £  assumed  small,  (ii)  the  coefficient  of  y~  assumed  small,  but  in  each 
case  the  convergence  of  the  results  gave  rise  to  doubt  as  to  the  legitimacy  of  the 

*  This  integral  and  also  that  in  equation  (89)  are  obtained  from  equation  (86)  by  means  of  the 
integrating  factor  yn.      It  may  be  that  with  other  integrating  factors  other  integral  harmonics  might 
afford  integrable  forms,  but  these  have  yet  to  be  found. 

f  At  the  close  of  §  18  of  the  memoir  cited  above. 

rti 

*  It  is  the  absence  of  the  linear,  term  f  ^  from  the  equation  for  the  boundary  obtained  by  Saint- 
Venant  which  forms  the  basis  of  the  criticism  above. 
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assumptions  by  failing  to  give  appropriate  values  for  the  quantities  supposed  small. 
Moreover  the  restrictions  imposed  on  the  constants  of  equation  (90)  would  be  much 
more  stringent  than  those  affecting  Saint-  Venarit's  symmetrical  boundary  (87)  and, 
even  if  the  analytical  difficulties  could  be  overcome,  there  would  be  correspondingly 
less  interest  in  the  solution. 

§21.     A  general  method  of  solution  by  means  of  conjugate  functions. 
The  case  of  the  circular  sector  cross-section  with  which  this  paper  is  chiefly 
concerned  might  be  regarded  as  a  solution  in  terms  of  the  conjugate  functions  given 

by  the  relation 

a  +  i$  =  log  (x  +  iy), 

and  it  is  clear  that  the  use  of  other  sets  of  orthogonal  curvilinear  coordinates  given 
by  the  general  relation 


might  give  rise  to  solutions  which  would  come  within  the  processes  of  analysis  to  an 
extent  sufficient  to  permit  of  arithmetical  reduction. 

With  a  view  to  simplifying  the  analysis  we  may  take* 


(x,  y)  =  -  rxy  -)-        (  1  +  1>?)  (xf  -  ±xs)  +  (a  -  a")  y  -  (&.  - 


........................  (91) 

where  $  (x,  y)  satisfies 

d'2(f>    c°-(f> 

3^  +  ^  =  0. 

ox      oy- 

With  this  form  of  x  it  is  seen  that  the  boundary  equation 


takes  the  simpler  form 

y-dx  =  0  ............  (92). 


Now,  by  the  properties  of  conjugate  functions, 

d(f)      d<f)da 


Sx      dadx      dfldx      da  ox      9/3  dy 


d(f>da      d(f>d/3      d(f>da. 
=  = 


and  along  the  contour  given  by  a  =  constant 

da  ,       da, 

—  dx  +  —  dy  =  0. 

dx          Sy 

*  In  order  to  avoid  confusion  of  the  slide  constant  ft  of  the  earlier  part  of  this  memoir  with  the 
conjugate  function  symbol  /3  we  now  write  the  former  /3  as  /?.,. 

8—2 
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Hence  if  the  contour  a  =  constant  is  to  be  a  part  of  the  boundary  of  the  cross-section 
of  the  beam,  we  must  have,  from  equation  (92), 

da 

Sd>  Sx  f  W, 

a    =7TT^ 7TTJ  2r2/  +  T?  r  \\ 

ca.      (oa\       loa\  M/i 


along  that  part  of  the  boundary,  and  similarly  if  ft  =  constant  is  to  be  a  part  of  the 
boundary  go 


along  that  part  of  the  boundaiy. 

The  analytical  problem  thus  becomes  : 
To  find  a  function  <j>  (x,  y)  such  that 


and,  therefore,  equivalently  by  the  properties  of  conjugate  functions, 

dx1      dy'2 

throughout  the  region  enclosed  by  a  closed  curve  composed  of  parts  of  the  contour 
lines  given  by  taking  a  or  /3  to  be  constant,  and  such  that  ^-  and    Q  assume  the 

values  given  in  equations  (93  a)  and  (936)  over  the  parts  of  the  boundary  where 
a  =  constant  and  ft  =  constant  respectively. 

The  algebraical  difficulties  that  arise  are  1°,  to  find  the  inverse  expressions  for 
x  and  y  in  terms  of  a  and  ft  for  substitution  in  equations  (93 a)  and  (936)  and  2°,  the 
determination  of  the  Fourier  solution  for  <jt  which  will  satisfy  equations  (93  a)  and 
(936)  in  the  way  required.  A  failure  to  surmount  these  difficulties  analytically  would 
not  however  render  the  solution  impossible  for  we  can  always  fall  back  on  arithmetical 
methods  to  solve  any  particularly  interesting  case.  Indeed  when  it  comes  to  the 
application  of  the  Conditions  I  and  II  (Equations  (23)  and  (24)  of  §  4),  it  will  be 
necessary  in  all  but  the  simplest  cases  to  use  arithmetical  or  mechanical  methods  of 
quadrature  to  determine  the  value  of  the  torsion  constant  r. 

In  the  case  of  symmetry  perpendicular  to  the  plane  of  flexure  Condition  I  is 
identically  satisfied  and  Condition  II  becomes : 

r  X  f  tifdxdf/  =  -        |~  I  f  (r)y>  -  ( I  +  17)  sfy)  dxdy- 2(1+1))  ly] 
J  J  z-/jj  U  J  J 

-  £  I  <f>  (x  cos  (yv)  —  y  cos  (xv})  ds (94). 

The  integration  with  regard  to  ds  must  then  be  expressed  as  an  integration  with 
regard  to  a  or  ft  as  the  case  may  require. 
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§  22.     Particular  cases  of  the  solution  by  conjugate  functions. 
(i)     The  Lemniscate  functions*. 
These  are  given  by 


and  consist  of  the  families 


and 


f  I 

a  +  iB  =  log —, -,  +lop-     — TT — 

0  x  + 1  (y  +  c)  x  +  %  (y  —  c) 


._,  • 

tan  p 


•(95), 

.(96) 
.(97), 


the  latter  being  a  system  of  hyperbolas,  each  of  which  passes  through  the  limiting 
points  (0,  ±c)  or  foci  of  the  system  of  lemniscates  given  by  equation  (96). 


I  =  centroid  of  lemniscate. 


€„  =  centroid  of  sector.        F=  focus  of  lemniscate  (OF=  c,  OD  =  a  =  c 
FIG.  11. 


*  Compare  Lame,  Lefons  sur  les  coordonnfas  curvilignes,  Paris,  1859,  p.  217  et  seq.,  where  will  be 
found  (p.  223)  a  diagram  of  the  two  orthogonal  families  of  curves.  The  notation  used  here  differs  to  some 
extent  from  that  used  by  Lame ;  it  is  arranged  to  provide  a  continuous  system  of  /3-values  for  the  Fourier 
expansion.  The  positive  half  of  the  axis  of  x  and  origin  to  focus  of  the  axis  of  y  correspond  to  /3  =  0,  the 
axis  of  y  beyond  the  focus  is  /3  =  TT,  the  negative  portion  of  the  axis  of  x  and  the  focus  to  the  origin  of  the 
axis  of  y  are  both  /?  =  2ir.  The  lower  branches  of  the  hyperbolae  through  the  focus  are  provided  by  values 
of  |3  between  0  and  TT  and  the  upper  by  values  of  /3  =  TT  to  2-rr,  the  two  branches  meeting  at  the  focus  corre- 
spond to  if/  and  TT  +  <{/,  see  the  diagram  above. 
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When  a  =  0  equation  (96)  reduces  to 


+  2c2(ar-?/)  =  0  ...........................  (98). 

This  is  the  Cartesian  equation  of  the  ordinary  lemniscate  r2  =  2c2cos  20,  a  curve 
which  takes  the  form  of  a  figure  of  eight,  the  tangents  at  the  double  point  (0,  0)  being 
the  octant  lines  x  =  ±  y. 

We  propose  to  find  the  solution  of  the  flexure  of  a  beam  whose  cross-section  is 
one  loop  of  the  general  lemniscate,  the  direction  of  flexure  being  in  the  direction  of 
the  axis  of  a;  and  therefore  at  right  angles  to  the  line  of  symmetry  of  the  cross-section. 
It  is  easy  to  show  that  along  the  boundary  of  this  single  loop  we  must  take  ft  to  vary 
from  0  to  2Tr  and  thus  complete  a  circuit  of  the  loop. 


Let  v=l—  e  "cos ft,  u  =  -J\—2e~ 

then  «7c2  =  |  (u-v),  y*/c2=i(u  +  v) (99). 

d<h  x  (x*  +  rf  +  c2)  f             W 

and                           /  =  —  -  ;r7-2yr-ir~  {2ry  ' 

//.-,  O   /  '**2  _L.  «J»I                      ^ 

'id.  ^/ 

which  leads  to 


We3  — 

1) (100), 


u0  and  va  being  the  values  of  u  and  v  when  a  =  a0,  the  value  of  a  for  the  required  loop. 
This  is  the  boundary  condition. 

Condition  I  is  identically  satisfied.     Condition  II  leads  to 

((  u.W  f  C 

2p.r  I    tfdxdy  =  Zjjfj  I  I  (rjf  -  (1  +  17)  a*y)  dxdy 

—  Wy  —  /u,  |<^a=a0(«cos  (yv)  —  y  cos  (xv))  ds  ......  (lOla); 

see  (24),  p.  12,  and  below  (37),  p.  17. 
If  F=  (x-2  +  f  -  c2)2  +  4c2cc2  -  c4e-2"  : 

dF/dx  =  4x  (x2  +  y1-  +c2),     dF/dy  =  4y  (cc2  +  f  -  c2), 

cos  (Xv]  =  dF/dx/((dF/dxY 

•    -i  /     \     y 

and  similarly  cos  (yv)  = 

•-     yce 

2xu  c  <u2  —          , 

Thus         {xcos  (yv)  —  ycos  (xv}\  ds=  —-r=   =^eads—  ---  ,  —      eds. 

*       > 


i  i  \«     /  7  \2 
But         (ds)-  =  (dxf 


4  w  u —  v  4  ?< 

or,  ds  =  lj  —r-   dft. 
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Accordingly  </>a  =  a0  (x  cos  (yv}  —  y  cos  (xv))  ds=  —  -       <£a=c 

Jo  ^  Jo 

=  ~  ^      <£«=<>« 


u 
sin, 


-  2e'ao  cos  £  +  e'2"" 
Thus  (tOla)  becomes 

W 


rxii 

_rf^ (1016). 

-2*1  ' 


Now  a  solution  of  -j*  +  --i  _ 

- 


is  <f>  =  '£  (An  sinh  na  +  Bn  cosh  na)  sin  n(ft--rr)  ............  (102  a), 

Inside  the  loop  a  =  0,  a  ranges  from  0  to  +  oo  at  the  focus  (0,  c),  but  w  and 
therefore  <^>  must  be  finite  ;  it  follows  therefore  that  An  +  Bn  must  be  zero  and  our 
solution  will  take  the  form 

3-Tr))  ........................  (1026), 


-'M°smn(p-Tr))  ............  (103). 

We  have  accordingly,  from  (100)  and  (103), 

LT//V- 

S  (nAne-»°°siu  n(/3  -  *))  =  ^f(/3)  +  fijj  *  (ft), 
where      f(p)  =  7<  -  <  (1  +  w.-1),    ^  (/«)  =  («.  ~  ».-  2i?v0)  4T«»  Z  «.l  (1  +  «,-). 

If  ATC''  + 


.,>lan  /"-2ir 

it  follows  that  ^n'=  -      -T-       I   f(B)smn(B-Tr)d/3, 

mr  +  sm  mrJ0 

An"= 


rnr  +  sin  mr  J  0 

These  integrals  must  be  obtained  by  quadratures  and  the  values  of  An'  and  AH"  then 
inserted  in  the  )8-integral  term  of  the  torsion  equation  (1016),  where  again  for  the 
general  case  quadrature  is  needful.  The  proper  values  for  n  can  only  be  determined 
when  we  have  found  the  expressions  f(B)  and  i//  (B)  for  the  particular  case. 

We  will  illustrate  the  method  on  the  relatively  simple  problem  of  the  torsion 
flexure  of  a  beam,  the  section  of  which  is  the  loop  of  an  ordinary  lemniscate.  In 
this  case  an  =  0  and 


J 


7V   T=r-  =        =°  cos 

o  v  1  —  2e  ~  "o  cos  B  +  e  '  2a«          Jo 

Now  in  this  case 

v0  =  2  sin2  J/3,     MO  =  2  sin  |/3, 
and  we  find 

x  =  csin^(ir-B)  N/2sm|/3,    y  =  c  cos  ^  (IT  -  B)  J'2  sin  |/3  ......  (104) 
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over  the  contour  a  =  0,  and  accordingly 


-y  - 
da/a=o 


=  —  I  TC2  (sin  /3  +  cos 


_ 
2  sin  |/8{cos  1(3/8  -  TT)  +  (£  +  17)  (sin  i  (£  -  TT)  -  sin  £  (5/8  -  TT))}  ......  (1046). 

Thus        /(/8)  =  sin/8  +  cosi/8, 

i/,  (£)  =  %/2smi0{cos  1  (3/8  -  TT)  +  (\  +  r?)  (sin  i  (/8-  TT)  -  sin  ^  (5/8  -  TT))}. 
A  suitable  form  of  <£  is 

<£  =  |TC'  (e  -  a  sin  /8  +  2e  ~  4"-  cos  £/8) 


e-'!asinn(/8-7r)    ......  (105). 

This  will  satisfy  (1046)  provided  we  take 

i  r2"  / 

An'"  =  -         s/2  sin  |/8  cos  £  (3/8  -  TT)  sin  n(ft-n)  d/3 
mr  J  0 

1        T2^       _ 

Aniv  =  -         s/2  sin  |/8  (sin  i  (0  -  IT)  -  sin  J(5/8  -  77))  sin  n  (/8  -  TT) 

?i7T  J  o 

These  may  be  shown  to  be 


...(106). 


A  '"  — 

A-n      —   — 


mr 


cos         sn    /    sn 


Aniv  =  —  r  (2  cos  |/8')*  sin  f/8'  sin  n/3'd 
mr  Jo 


.(107). 


These  integrals  were  found  by  quadratures,  using  one  of  Sheppard's  formulae  (Bio- 
metrika,  Vol.  I,  p.  276,  Case  (i)  (c)),  the  subject  of  integration  being  tabled  to  every  5°. 
After  considerable  labour  there  resulted  : 


n 

^n'" 

-*„'" 

«(-!)» 
•     4n2-l 

*n  =  Mn"'  +  Mn" 

"»'=;Mnlv 

1 

-1-125,0142 

+  1-500,0012 

-  -333,3333 

-  -093,7535 

+  -375,0003 

2 

+    -031,2632 

+    -187,4978 

+  -133,3333 

+  -031,2530 

+  -046,8745 

3 

•013,0284 

-    -020,8406 

-  -085,7143 

-  -005,8622 

-•005,2102 

4 

+    -006,8490 

+    -005,8517 

+  -063,4921 

+  -002,4437 

+  •001,4629 

5 

-    -004,0996 

-    -002,3637 

-  -050,5051 

-  -001,3204 

-  -000,5909 

6 

+    -002,6992 

+    -001,1255 

+  -041,9580 

+  -000,8155 

+  -000,2814 

7 

•001,8626 

-    -000,6602 

-  -035,8974 

-•000,5482          --000,1650 

Now      <f>a=o  =  ITC"  (sin  /8  +  2  cos  A/8)  +  2  (A,' 


We3 

:  sin  n  (/3  —  77) 


and 


a=o  cos  1/8^/8  =  ^c2  (f  +  27T) 


°°  r 

S     ( 

i     |^ 
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Equation  (1016)  for  the  torsion  now  becomes 


W 

lj)y-l+:>yxy-2l+ 

*)  ......  (108). 

\\y2dxdy  =  ±c4(\Tr  +  %},      \  \y3  dxdy  =  \%TTC\      \\tfydxdy  =  -faird>, 


But  for  the  common  lemniscate 


thus  we  have 


We 

T  (-  7853,9816)  =  ---{'002,4468  +  '121,3194^}, 


We 
T=  -Wr  {'003,1  154  +  -154,46877;}  .....................  (109) 


=  _  ^,"{-002,203  +  -109,22597/},  if  a  =  « 


the  breadth  of  the  lemniscate.     If  77  =  5, 


r=  -  Sf.x  -02951. 


Thus  if  the  lengthy  analysis  be  without  slip,  there  is  a  small  negative  torsion,  the  first 
we  have  come  across. 

We  can  now  complete  the  solution.     We  will  write 

We. 


and  we  have 

We3 


e-«-  sn  » 
sn  j   +  e-a  cos/e1  +  e2Ti)  -  S  {(/c,,  +  Kn'rj)e~  ""sin 


the  numerical  values  of  «„,  «,/  being  given  on   p.  64.      Hence  by  equations  (20) 
and  (91) 


Jxy 
El  Ll     £37?^  V  c*      ^        81U  P~e       cos 


n 

c  i 

and  y  will  be  as  in  equation  (20),  p.  12. 


sn  n       -  TT     - 
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We  have  now  to  determine  the  fixing  of  the  clamped  end  ;  let  us  suppose  the 
element  at  the  centroid  fixed.     Then,  by  p.  42,  as  =  P'  =  0  and 

W 


WIL  _2 

Ei^y 

w 

Hence  « (0,  y)  =  —  ryz  +  -~. 


W  nj,_,    3} 

.W 


gives  the  linear  term  in  the  deflection  as  soon  as  fts  is  found. 

We  have  now  to  find   ,    along  the  axis  of  symmetry.    Clearly  -^  =  0,  when  x  —  0, 


dft 

7—  =  --  2_       c      when  x  =  0. 
ax  c  —  y 

But  e-a  =  (c"  —  2/2)/c2  for  y  <  c,  =  (y*  —  e2)/ca  for  y>c,  and  at  the  centroid  the  former 
must  be  used  for  y  =  -^Trc.  Further  /3  =  0  for  x  —  0,  from  the  origin  to  the  focus,  and 
—  TT  beyond  the  focus,  i.e.  y  =  c.  Thus  between  origin  and  focus 

dw 


-     S  2n  (  -  1  )»  e  -  (»  -  1)  -  (KH  +  K/T,)    -  /8S. 
v  i 

Further,  by  equation  (20)  along  x  =  0, 


Hence  adding 

I  xz\  div     du      Wc'V,  .y        y" 

-     1  =  0=      ,       +   -y-  =      7-iT      (CJ+ C..T7)  - +77  — 

\  /ot  /y<c      dx      dz       El  Lv  '  c      '  G- 

...(110), 


which  enables  us  to  get  the  shearing  stress  from  the  origin  to  the  focus.  Beyond 
the  focus  /3  =  TT,  thus  cos/3  changes  sign,  and  c3  —  if=  -c2e-Q,  thus  the  term  |e-asiny8 
in  w  remains  the  same  as  before  in  dw/dx,  but  the  last  term  changes  and  we  have 


xz 
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At  the  focus  a  =  oo  and  the  summation  reduces  to  a  single  term  the  same  for  both ; 
thus  we  have 


F  /1t~t 


Wc\ 
EI\ 

We* 


1+2^)77] 

-1-904,469377] (U1) 

Wa* 
•29173,  if  77  =  ^,  =  ^y-x  -14586,  ifaW2c. 


At  the  centroid  we  have 


and 


,  ,  x -, 

^"          "  ~~ 


=  '3831,4973. 


But 


.rz 


VFic3  f, 

•r=°  =  ^V     V€' 
i/  =  p       -fc^   L 


2     ' 


*\ 


or,  the  curve  of  deflection  is 


as  it  should  be.     It  only  remains  to  find  /3,  and  s0. 
Numerically  we  have 


fj-/y=y 


.(112), 


71 

iTne-t"-1)"!-!)" 

/rnx47rrei?-(n-1)a(-l)" 

/rB'xJirne-lH-1la(-l)n 

1 

-1-570,7963 

+  •147,2677 

-•589,0491 

2 

+  1-203,7003 

+  •037,6192 

+  -056,4228 

3 

-    -691,7963 

+  -004,0554 

+  -003,6044 

4 

+    -353,4154 

+  -000,8636 

+  -000,5170 

5 

•    -169,2637 

+  -000,2235 

+  -000,1000 

6 

+    -077,8238 

+  -000,0635 

+  -000,0219 

7 

-    -036,3589 

+  -000,0199 

+  -000,0060 

These  values  lead  us  to 
Wo\ 


, 425177- -190,1128 +  -528,3770^] 


We*, 


=  «[ -'190,1 128  +  -83'5,8021 77]  ... 

Wo? 

=  -=,  X  '009,544, 
Jiil 

if  77  =  1  and  a=«/2c  be  the  breadth  of  the  lemniscate  loop. 


.(113) 


Further, 


We 


sa=      ^[1-266,546777  --187,6660] 


(114) 


Wn" 
=    —  x  -06449, 


if  77  =  ^  and  o=v2c'  be  tlie  breadth  of  the  lemniscate. 
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Thus  the  deflection  will  be 


The  linear  term  is  only  about  ^-th  of  the  linear  term  in  a  circular  sector  of  angle  90°. 
This  is  due  to  the  change  in  the  torsion  from  a  sensible  positive  to  a  slight  negative 
value. 

It  remains  to  investigate  the  vertical  shear  along  the  central  radius.  In  working 
this  out  attention  must  be  paid  to  the  change  in  form  in  the  shear  after  passing  the 
focus.  The  following  table  gives  the  chief  results  : 

Shears  parallel  to  load  along  Axis  of  Symmetry. 


ylc. 

yla 

£/** 

1=4 

-.    /  lI//*2 

"/IT 

xz/(mean  shear)  -  1 

0 

0 

0 

0 

-  1 

•1 

•0707 

-  -019,541  +    -038,2897; 

-  -0100 

-1-135 

•2 

•1414 

-•037,932+    -087,57377 

-  -0160 

-1-216 

•4 

•2828 

-•067,834+    -222,5627; 

-•0122 

-  1-164 

•6 

•4243 

-•086,317+    -410,126); 

+  •0162 

-    -782 

4*" 

•5554 

-  -093,833  +    -633,2737; 

+  -0645 

•131 

1 

•7071 

-•092,196+    -952,235); 

+  -H59 

+    -966 

1-15 

•8132 

-•084,637  +  1-214,414)7 

+  -2190 

+  1-951 

1-30 

•9192 

-•07  1,986  +  1-509,  133); 

+  -3053 

+  3-114 

V2 

1-0 

-•060,879  +  1-754,88477 

+  -3778 

+  4-091 

The  last  column  is  plotted  on  the  diagram,  p.  61.  It  is  clear  that  from  the  origin 
onwards  we  have  a  small  negative  shear  and  ultimately  a  moderately  large  positive 
shear,  reaching  a  maximum  at  the  butt  end  of  the  section.  The  form  of  the  dis- 
tribution differs  widely  from  that  of  the  90°  circular  sector.  We  have  not  succeeded 
in  discovering  any  physical  reason  for  the  torsion  being  negative  in  the  case  of  the 
lemniscate  and  positive  in  the  case  of  the  circular  sector.  On  the  other  hand,  we 
have  not  discovered  after  much  investigation  any  slip  in  the  analysis.  We  publish  the 
investigation  in  the  hope  that  a  fresher  eye  may  discover  our  slip  or  verify  the  analysis, 
or  again  find  other  sections  with  negative  torsion.  Should  such  sections  actually 
exist,  the  quick  transition  in  sign  with  changes  of  not  very  marked  character  in  the 
asymmetry  might  suggest  that  the  edge  of  attack  in  a  propeller-blade  may  be  of  a 
rather  unstable  nature.  Possibly  when  time  is  available  for  the  sections  suggested 
below  to  be  worked  out  numerically,  light  will  be  thrown  on  this  interesting 
question. 

(ii)     The  coaxial  circle  functions*. 

The  equation  x  +  i  (y  +  c) 

a-Mp  =  log         .;  (116) 

x  +  i  (i/  —  c) 

*  '  Systems  cylindrique  bicirculaire'  of  Ltmn'; :  see  his  Le.fons  sur  leg  coordonnee.s  curmliynes,  p.  199, 
Paris,  1859. 
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gives  rise  to  the  two  families  of  circles 

y?  +  (y  —  c  coth  a)2  =  c3  cosech2  a  .....................  (117), 

(x  —  c  cot  B)'*  +  y"  =  c'  cosec2  (3    ................  .....  (118), 

the  a-family  being  a  system  of  coaxial  circles  having  real  limiting  points  at  (0,  ±c) 
and  the  axis  of  x  as  radical  axis,  and  the  /8-family  being  the  orthogonal  system  of 
coaxial  circles. 

The  interesting  case  from  the  present  point  of  view  is  the  boundary  given  by 
taking  two  circles  of  the  a-family,  one  being  wholly  inside  the  other*.  The  beam 
then  is  a  hollow  shaft  with  an  eccentric  cavity,  and  thus  possesses  the  requisite 
uni-axial  symmetry. 

The  expressions  we  require  for  substitution  are  : 


csin  B 
x  = 


cosh  a  —  cos  B 

c  sinh  a 

11  = — 

9     cosh  a  —  cos  B 


•(119), 


dx 


+ 


da 

sinh  a  sin  B    -j 

dx 
do.  _ 

c 
1  —  cosh  a  cos  B 

ty 

*\  \  c, 
day 

c 
(cosh  a  —  cos  By 

(120), 


and  it  is  easily  seen  that  a  complete  circuit  of  any  circle  of  the  a-family  is  given  by 
taking  B  to  vary  from  —IT  to  +  TT. 

Equation  (93)  becomes  for  this  case 

We?  (I  +  77)  sinh  asin3/3  —  77  sinh3asin/3 
(cosh  a  —  cos  B)3 

..................  (121), 


d(f>  _  ,     sinh2asin/3 

3a  *  (cosh  a  —  cos  3 


El 


and  the  form  of  the  solution  will  bef 
<        2 


smh  n  (a,  —  a2) 

where  the  constants  An  and  Bn  are  determined  by  Fourier's  method  from  the  equations 
obtained  by  substituting  at  and  a2  for  a  in  equation  (121). 

Condition  I  is  again  identically  satisfied.     Condition  II  becomes 

w 


r   sin  £  cosh  a  ...... 

J  _„  cosh  a  —  cosBT 

*  This  problem  is  suggested  as  a  soluble  one  by  Love  (Elasticity,  Second  Edition,  §  232  (e),  p.  325). 
The  corresponding  torsion  problem  has  been  investigated  by  H.  M.  Macdonald,  Proc.  Camb.  Phil.  Soc., 
Vol.  viii.,  p.  62,  1893. 

t  Lame,  loc.  cit.,  p.  184. 
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In  the  /3-integral  a  in  <f>  and  cosh  a  must  be  given  the  two  selected  values  and  the 
integrations  round  the  two  circles  taken  in  opposite  senses.  As  in  the  case  of  the 
lemniscate  <f>  will  contain  series  terms  with  both  r  and  WjEI  as  coefficients. 

To  get  the  solution  for  a  beam  resembling  in  cross-section  the  blade  of  a  propeller 
we  might  use  the  curvilinear  quadrilateral  enclosed  between  the  circular  arcs  given  by 
dj,  a.2,  y8,,  —  /61;  but  in  the  sector  of  an  annulus  we  have  such  a  similar  figure  that 
there  would  be  little  but  the  pleasure  of  confirmation  in  following  out  this  case. 

(iii)     The  confocal  conic  function**  '. 

The  appropriate  systems  for  our  problem  are  given  by 


)   ........................  (124), 

and  consist  of  the  family  of  ellipses  with  foci  at  (0,  ±  c) 

-y'2  i/2 

+    _if  _  =  c*  ..  ..(125) 

sinh2  a      cosh2  a 

and  the  family  of  hyperbolas  with  the  same  foci 

*  if 

+      7    =  C'">" 


The  Cartesian  coordinates  of  the  system  are 

acos/3 


...........................  (127) 

asin/SJ 

• 

and  it  is  seen  that  the  part  of  the  axis  of  y  lying  between  the  foci  may  be  taken  as 
the  ellipse  a  =  0,  and  the  part  lying  outside  of  the  foci  as  the  hyperbola  /3  =  ^7r,  and 
the  positive  axis  of  x  as  the  hyperbola  /3  =  0  and  the  negative  axis  of  x  as  the  hyper- 
bola {$  —  IT. 

From  the  point  of  view  of  the  flexure  problem,  and  especially  in  relation  to  the 
problem  of  the  propeller-blade,  it  would  appear  that  the  most  interesting  case  would 
be  the  boundary  enclosed  between  an  ellipse  a0  and  one  limb  of  a  hyperbola  given 

by  taking  $  =  --—$,,  where  /?„  has  an  appropriately  small  value  so  that  the  eccen- 

Zt 

tricity  of  the  hyperbola  will  be  fairly  large.  As  in  all  cases,  however,  where  composite 
boundaries  are  made  up  from  systems  of  orthotomic  curves,  there  is  the  disadvantage 
of  the  right-angled  corners  which  do  not  appear  in  the  section  of  a  propeller-blade. 

For  this  case,  where  —  a,  <  a  <  a0  and  -  —  /30</3<-  +  /60,  the  boundary  conditions 

Zt  £ 

derived  from  equation  (93)  are  easily  found  to  be 

3<f>  W<? 

r—  =  2TC2  cosh2  a  sin  /6  cos  ft  +  -p-j-  ((l+if)  sinh2  a  cosh  a  cos3  ft  —  rj  cosh3  a  sin2  /3  cos 


;-  ........................  (128) 

*  'Cylindres  homofocaux'  of  Lam£,  loc.  cit.,  p.  195. 
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and 

~  =  —  2rc2  sinh  a  cosh  a  sin2  B  — ^7-  (( 1  +  >?)  sinh3  a  sin  /3  cos2 18  —  17  sinh  a  cosh'2  a  sin:i  /3) 
dp  tLl 

(129), 

and  the  solution  for  <£  will  be  of  the  form 

.   .   mr  I  n     TT     0\        nir  ,  \      n     •  -Ln^  /  n>a  I  a     ™  ,  0\ 

$  =  An  sinh  -    [  /3  -  -  +  j80   cos  —  (a  -  a,)  +  Bn  smh  -=-  (a  -  a0)  cos  -=-  h8  -  -  +  0. ; 

«o    \          ^  /  ao  Po  ft  \      :    *  / 

.     (130), 

where  the  constants  .4n  are  derived  by  Fourier's  method  from  equation  (129)  for 

B  =  -±Bs  and  the  constants  jB,,  from  equation  (128)  for  a  =  a0. 

Lt 


an 


Turning  to  (94)  the   \<$>(x  cos  (yv)  —  y  cos  (*v))  ds=  —  ^c2  I  sin  2/8^>o0<^/8  along 

a  =  a0  boundary,  and  to  +^c2  I  sinh  2a<f>p1da  along  a  p\  boundary.     Thus  the  equation 
for  the  torsion  is  provided  as  soon  as  An  and  Bn  have  been  determined. 

It  is  evident  that  this  solution  can  be  completely  set  forth  ready  for  arithmetical 
calculation  without  methods  of  mechanical  integration. 

(iv)     The  rosette  functions. 

We  shall  lastly  consider  the  system  of  conjugate  families  given  by 

a  +  ift=    (al^  ..(131). 

(x  +  iy)m 

When  polar  coordinates  (a5  =  rsin#,  y  =  rcos0)*  are  used  the  equations  of  the 


two  families  are  found  to  be 


*  /» 

a  =  —  -  cos  mv 
r 


and  we  easily  deduce  that 

-  —        /I  ft 

x  =  a  (cr  +  /32)   2'"  sin  (  -  tan  1 
\m  a 

11  =  a  (a2  +  B*Y  *'"  cos  f  —  tan"1  - 

\  "?y)  rt 

and  that 

9a 

9x  yiH-i 

..  .    '    /0  .,  = sin  (m  +  1 )  6 

daV  .  fda\-         mam 

(fa    ' 


m  m  a/  \m 


*  An  unusual  form  of  polar  coordinates  has  been  adopted  in  order  to  preserve  uniformity  with  the 
notation  given  in  the  earlier  portion  of  the  paper. 
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The  typical  curve  of  the  a-family  is  a  rosette  composed  of  petals  radiating  from 

_j\_ 
the  origin  and  enclosed  within  a  circle  of  radius  aa  m,  which  is  touched  by  each  petal. 

Let  us  take  the  rosette  a  =  1  and  let  us  consider  the  single  petal  which  results  from 

the  principal  values  of  tan^/S  being  used,  that  is  for  —  -  ^tan"1  /3<  +  -.     This  petal 

z  z 


is  a  closed  curve  symmetrical  about  the  axis  of  y  and  having  a  node  of  angle  —  at  the 

origin.     With  a  suitable  choice  of  m  we  thus  have  a  curve  which  would  resemble  fairly 
closely  the  cross-section  of  a  propeller-blade,  and  it  is  with  this  in  view  that  we  proceed 
to  give  a  formal  solution  of  the  flexure  problem  of  a  beam  having  such  a  cross-section. 
The  boundary  equation  (93)  becomes  after  reduction 

1    ii  t 

m  cos   —  tan 


-  —  -^Y 

3a/a=1          L  m  \m  I       El       m 

sin2  ^tan-'/3)^ 

........................  (135), 

and  this  is  to  be  satisfied  for  the  whole  of  the  range  —  <x>  <  /3  <  oo  ,  a  range  which 
renders  a  Fourier  series  solution  inadequate. 
Consider  the  integral 


1   f°°     fx  fe-«      e~w(aL'l) 

=  -  &(*){•  -cosu(v-0))dudv     ...(136). 

7Tj-x  Jo  \   U  U  'I 


The  integrand  is  everywhere  finite  if  F  (v)  is  finite  for  —  oo  <  v  <  oo  ,  and  it  will 
be  assumed  for  the  present  that  the  integral  is  convergent*  for  the  range  0  $  a  $  1, 
-  oo  <  /?  <  oo  .     Also  it  satisfies  the  differential  equation 


and  we  have 

/3r/>\  1     /"X  fx 

(r-M      =  -l      F(v)  I    cos  u  (v  —  B)dudv (137) 

\3a/a=1      »•/_<,         'Jo 

=  F  (/3),  by  Fourier's  integral  theorem. 

If,  therefore,  for  F  (v)  in  equation  (136)  we  substitute  the  expression  for 

(a,  V)\ 


\      d<*-       /a=i 
*  It  does  not  appear  to  have  been  noticed  by  Fourier  that  the  integral  I  —       -  -  cos  u  (v  -  /?)  du  is 

divergent  for  any  range  of  integration  which  includes  the  point  u  —  0.  Compare  the  Theorie  analytique 
de  la  Chaleur,  §  421,  Darboux's  edition  of  the  (Euvres  de  Fourier,  Vol.  I.,  p.  507,  Paris,  1888,  where  the 
integral 

/°°  f"  dp 

daf(a)  I      —  cos  (px  —  pa)  (em  +  e~pv) 
-00  J -vt    P 

is  considered  without  any  mention  of  convergence ;  this  has  probably  been  often  remarked  before.  The 
convergency  of  the  integral  (123)  above  is  obtained  by  the  insertion  of  the  constant  I  I  F(v)—-dudt> 

J  -a,  JO  U 

which  is  also  infinite. 
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from  equation  (135)  we  obtain 


0     \         0a  «=1 

as  the  formal  solution  of  our  flexure  problem,  \vhen  combined  with  equation  (91). 
The  expression  for  (j)  (a,  ft)  takes  a  simpler  form  when  we  put  v  =  tan  w,  ft  =  tan  i/». 
The  value  of 


I  {x  cos  (yv)  —  y  cos  (xv)}  (f>ds 

in  the  torsion  equation  (94)  reduces  to  the  simple  expression 

1     ,  f+*    (<k-i)/3     ,0         «2  f+x,  1 

n«2  ^^=-2         <l>^d~        — 

"(i+^v5"  (i+/?r 

--*-|     '  I      r*F(v)s™u(v- 

Z7J-J   _x  J  o    J-ao 

which  completes  the  formal  solution,  but  is  an  expression  hardly  likely  to  prove  of  vahie. 
The  portion  of  the  integral  (136)  which  concerns  u  depends  essentially  on  the 
function  log  ((a  —  l)'J  +  (v  —  ft)'),   and  it  is  suggested   that  a   variant  formal  solution 
would  be  obtained  by  taking 

(139). 


Assuming   convergence,   we    can    easily    show    that    this   integral   satisfies  the 
differential  equation 

and  that 

O      =F(®    (140)> 

so  that  as  before  we  can  take  as  solution 


The  question  of  the  convergence  of  the  integrals  Used  in  this  discussion  clearly 
requires  further  investigation.  At  first  sight  it  might  appear  that  the  integral  in 
(139)  is  divergent,  seeing  that  logarithmic  term  contains  (v  —  ftf,  but  it  is  to  be 
remembered  that  F  (v)  is  in  our  case  an  odd  function  of  v,  so  that  it  is  not  decided 
that 

f      *»  1  °g  {(«  -  1  )2  +  (v  ~  PY\  dv  =  Lt    f  "   F  (v)  log  {(a  -  1 Y  +  (v  -  ft)*}  dv 

J   — 00  V-*-00    J    — V 

is  infinite.     There  is  further  to  be  decided  the  uniform  convergence  of  the  integrals 
when  ft  =  ±  oo  .     These  points  are  left  over  for  future  discussion. 

The  crude  formal  results  given  in  the  latter  part  of  this  section  may  prove  on 
further  investigation  to  be  unworkable,  or  at  least  unprofitable,  and  are  published 
here  only  by  way  of  indicating  the  lines  along  which  it  is  intended  to  continue  the 
work  when  opportunity  offers. 
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CONCLUSION. 

To  the  engineer  and  especially  to  the  designer  of  aeroplanes  the  chief  interest 
of  this  work  will  lie  in  its  application  to  the  problem  of  the  change  of  the  '  angle  of 
attack  '  that  arises  in  the  blades  of  a  propeller  rotating  at  high  speed,  and  it  may  be 
well  to  state  here  briefly  the  results  that  seern  to  affect  that  problem.  It  is  first  of 
all  to  be  emphasized  that  the  solutions  given  in  the  memoir  are  solutions  of  problems 
of  much  less  generality  than  those  of  an  actual  propeller-blade  and  that  the  results 
can  therefore  be  regarded  only  as  first  approximations  to  what  is  wanted.  Thus  the 
beams  considered  are  uniform  in  cross-section  throughout  their  length  and  that  cross- 
section  is  in  every  case  symmetrical  about  an  axis  which  is  perpendicular  to  the  plane 
of  flexure.  In  beams  of  cross-sections  having  no  axis  of  symmetry  the  torsion  effect 
would  presumably  be  of  greater  importance*.  Moreover  the  loading  is  not  continuous 
as  it  is  in  an  actual  propeller  but  consists  of  a  single  load  acting  at  the  end  of  the 
beam.  With  these  reservations  we  can  state  the  result  : 

For  a  beam  of  isotropic  material  and  of  uniform  uni-symmetrical  cross-section  of 
shape  resembling  the  cross-section  of  a  propeller-blade^  flexed  by  a  load  of  weight 
W  applied  at  one  end  in  a  direction  perpendicular  to  the  axis  of  symmetry,  the 
torsion  per  unit  length  due  to  the  unbalanced  shear-couple  may  be  taken  as  approxi- 

Wa 
mately  t  equal  to  '28  -^y  ,  when  E  is  the  Young's  modulus  of  the  material,  a  is  the 


breadth  of  the  beam  and  /  is  the  second  moment  of  the  cross-section  about  the  axis 
of  symmetry.  The  direction  of  the  torsion  is  such  that  the  thicker  side  of  the  beam 
twists  upwards  with  a  downwards  flexing  force. 

In  beams  of  cross-section  resembling  the  cross-section  of  a  propeller-blade  the 
torsion  due  to  anticlastic  curvature  is  in  the  opposite  sense  to  the  shear  torsion  and 
roughly  about  half  its  magnitude. 

If  the  material  of  the  beam  be  aeolotropic  we  have  seen  that  the  torsion  due  to 
the  unbalanced  shear-couple  may  be  increased  from  twice  to  six  times  the  isotropic 
value.  This  increase  is  largely  determined  by  the  ratio  E/'p.,  where  E  is  the  longi- 
tudinal stretch-modulus  and  p.  is  the  transverse  slide-modulus  parallel  to  the  direction 
of  loading,  but  it  is  far  from  being  directly  proportional  to  this  ratio. 

When,  however,  wood  is  the  material,  each  particular  beam  would  appear  to  have 
its  own  amount  of  torsion,  the  asymmetry  of  the  elastic  properties  of  the  material 
being  apparently  more  effective  than  the  asymmetry  of  shape.  To  the  designer  of 
aeroplane  propellers  this  should  be  of  great  interest  for,  as  mentioned  above  J  after 

*  The  amount  of  torsion  will,  of  course,  be  a  maximum  when  the  load  is  perpendicular  to  the  axis  of 
symmetry.     If  it  make  an  angle  if/  with  it,  only  the  component  W  sin  <j/  produces  torsion. 
t  For  the  calculations  on  which  this  estimate  is  basfd  see  §§10  and  13. 
t  In  §19. 


ON  THE  TORSION  PRODUCED  BY  FLEXURE  75 

the  description  of  an  experiment  which  brought  out  the  importance  of  this  asymmetry 
of  texture,  propellers  might  be  constructed  of  layers  of  wood  selected  and  arranged 
in  such  a  way  as  to  minimise  the  torsion  which  arises  from  flexure. 

Apart  from  the  pi'opeller-blade  approximation  the  most  interesting  case  theoreti- 
cally is  that  of  the  split  pipe  in  which  the  torsion  effect  is  so  pronounced.  We  do  not 
know  of  any  instance  in  practical  engineering  in  which  this  is  of  much  importance 
although  such  may  well  exist,  but  it  is  conceivable  that  some  practical  value  could 
be  found  for  the  phenomenon  in  such  a  pipe,  for  example  as  an  economical  form  of 
cantilever  weighing  machine  for  heavy  loads  or  in  circumstances  in  which  the  ordinary 
methods  would  be  inconvenient.  It  might  also  be  used  for  light  loads  as  a  letter- 
weighing  physical  toy.  Such  value  arises  from  the  fact  that  the  angle  of  torsion  is 
always  proportional  to  the  suspended  load  and  is  very  large  if  the  wall  of  the  pipe 
be  thin. 
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